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The  temperature  dependence  of  transport  properties  of  a  moderately 
dense  square  well  gas  is  studied  in  order  to  understand  the  effects  of 
attractive  forces  (particularly  bound  states).   The  quantum  cluster 
expansions  of  the  Oreen-Kubo  time  correlation  functions  for  the  thermal 
conductivity,  shear  viscosity,  and  self-diffusion  coefficients  are  given, 
and  exact  expressions  to  zeroth  (Roltzmann  level)  and  first  order  in  the 
density  are  obtained.   Specializing  to  Boltzmann  statistics  and  the 
classical  square  well  potential  allows  calculations  of  the  kinetic- 
potential  parts  of  the  first  densitv  correction.   Oood  agreement  with 
molecular  dynamics  results  is  found,  but  quantitative  differences  from 
real  fluids  are  observed.   Possible  reasons  for  the  discrepancies  are 
discussed.   The  important  contributions  to  the  remaining;  triple  collision 
parts  are  discussed.   A  brief  description  of  the  utility  and  limitations 
of  the  hard  sphere  model  is  eiven  for  comparison.   In  particular,  the 
dynamic  structure  factor  is  calculated  for  a  dense  fluid  of  hard  spheres 
and  compared  with  recent  neutron  scattering  data  for  Krypton  and  molecular 
dynamics  simulations. 


CHAPTER  I 
INTRODUCTION 

The  evaluation  of  the  equilibrium  anH  dynamic  properties  of  dense 

gases  and  liquids  relies  to  a  lartre  extent  on  theories  with  various 

degrees  of  phenomenology  or  uncontrolled  approximations ,  as  might  be 

expected  in  a  system  with  such  complicated  many-hodv  interactions.   For 

equilibrium  properties  (e.g.,  the  equation  of  state  and  other 

thermodynamic  functions)  there  exists  a  rigorous  theorv,  the  cluster 

expansion  method   developed  by  Ursell  and  Mayer,  for  finding  these 

properties  which  can  he  used  to  test  the  phenomenological  novels.   The 

cluster  expansion  generates  a  power  series  in  the  density  whose  first  few 

coefficients  are  in  a  form  suitable  for  practical  calculation.   Thus 

expressing  the  result  of  a  phenonenological  model  in  a  similar  series 

allows  one  to  determine  which  effects  the  model  properly  predicts  and 

those  for  which  it  falls.   For  dvnamical  properties  (e.g.,  transport 

coefficients)  there  is  at  present  no  such  theory  to  use  as  a  benchmark 

beyond  the  low  density  Holtzmarm  limit.   Although  the  attempts  of 

2 
Bogolubov  and  others   in  finding  a  systematic  derivation  of  the  Roltzmann 

transport  enuation  at  first  indicated  that  a  density  expansion  similar  to 

the  virial  expansion  might  he  found,  further  analvses  indicated  two 

complications.   Dynamic  properties  are  in  general  non-analytic  functions 

3 
of  the  density",  and  the  cluster  expansion  analogous  to  that  used  for 

equilibrium  properties  must  be  resumed  partially  to  five    the  correct 

expression  for  a  dynamical  property  as  a  function  of  the  density.   As  a 


consequence  there  has  been  no  rigorous  calculation  of  transport  properties 
beyonr!  the  low  density  boltzmann  result,  except  for  Senders  et  al.  , 
calculation  of  the  first  order  density  corrections  for  hard  sphere 
transport  coefficients.   The  main  portion  of  this  work  (Chanters  TIT,  IV, 
and  V)  will  concern  an  extension  of  Senders'  work  to  a  more  realistic 
potential,  so  t^at  the  effects  of  the  attractive  part  of  the  interaction 
on  the  transport  properties  may  be  examined.   It  is  known  that  the 
temperature  dependence  of  properties  for  both  real  fluids  and  potential 
models  with  attractive  parts  differ  strongly  from  hard  spheres 
particularly  at  low  temperatures.    However,  there  is  no  adequate  theory 
of  these  differences  outside  the  Boltzmann  limit.   For  exannle,  a  recent 
comparison  between  computer  simulation  and  phenomenological  models  show 
that  it  is  difficult  to  predict  even  the  correct  siffn  of  the  first  order 
correction  to  the  Boltzmann  transport  coefficients.    The  simplest 
interaction  which  incorporates  attractive  forces  is  the  square  well 
potential, 


<t>(r)  =  <  -e 


r    <   a 

a  <   r  <  T?a 


(1.1) 


0    r  >  Ra, 


where  r  is  the  relative  separation  of  the  particles  and  a   and  Ra  are  the 
inner  and  outer  radii  of  the  well.   The  transport  properties  to  be 
considered  are  the  shear  viscosity  (p),  thermal  conductivity  (X),  and 
self-diffusion  (nP  where  n  is  the  density).   [The  fourth  coefficient,  bulk 
viscositv  (k),  is  known  to  be  of  higher  order  in  the  density  than  the 
level  considered  here.]   The  square  well  potential  is  still  highly 


idealized  when  compared  to  a  more  realistic  potential  such  as  the  Lennard- 
Jones  model;  however,  the  square  well  model  retains  many  of  the 
simplifications  of  the  hard  sphere  model  that  allow  a  more  detailed 
treatment  of  the  calculations.   At  thp  outset  the  square  well  appears  to 
he  a  good  compromise  between  ease  of  calculation  and  introduction  of  the 
qualitative  features  associated  with  realistic  potential  models  (e.g., 
finite  collision  times  and  bound  states),  since  the  square  well  parameters 
may  he  adiusted  to  give  a  reasonably  stood  fit  to  the  second  virial 
coefficient  and  coefficient  of  shear  viscosity  of  the  Lennard-Jones  model 
(see  Figures  1.1  and  1.2).   However,  it  should  he  noted  that  the 
parameters  chosen  for  the  virial  coefficient  and  those  for  the  transport 
coefficients  are  in  general  quite  different.   Further  comment  on  the 
quality  of  the  square  well  potential  as  a  model  for  fluid  behavior  will  be 
given  in  Chapter  VI. 

Calculations  with  the  hard  sphere  model  can  generally  be  taken  to  a 
much  greater  depth  than  with  other  potentials.   For  example,  the  first 

non-analytic  density  corrections  to  Sergers'  results  h.^ve  been 

7  « 

calculated.    More  generally,  the  phenomenolotH  cal  Fnskog  model'  provides 

a  satisfactory  estimate  of  the  transport  coefficients  for  hard  spheres 

over  a  wide  ranre  of  densities.   Various  ad    hoc  modifications  of  the 

Q 

Fnskog  model  to  represent  other  potentials  have  also  been  suggested. 
Thus  hard  spheres  and  models  based  on  hard  spheres  are  still  used  in  many 
problems.   The  ability  of  hard  sphere  models  to  predict  the  behavior  of 
dense  gases  and  liquids  will  he  examined  in  Chapter  II.   There  it  will  be 
seen  that  agreement  of  hard  sphere  kinetic  theories  with  computer 
simulations  of  hard  spheres  is  generally  quite  good  over  a  range  of 


D 

r-i      t= 

TO   CM 


Figure  1.2.   Shear  viscosity,  n  ,  for  the  square  well  potential 

( )  and  the  Lennard-Jones  potential  (o  o  o)  as  a  function  of 

T*  =  kT/e  .   R  =  1.96;  eLJ/e0  =  -717;  a   Ja   =  1.143. 
Coefficients  are  normalized  to  hard  spheres. 
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densities,  but  that  quantitative  and  qualitative  differences  can  be  found 
when  compare^  to  experimental  data  or  otber  Dotentia]  models. 

The  transport  coefficients  will  be  evaluated  t^rou^h  their  Green-Kubc 

Q 

formulas  ;  anv  general  transport  coefficient,  £  =  nD,n,A,K,  can  be 
expressed  as  the  Laplace  transform  of  a  time-correlation  function: 


C  =  lim  -f  j   dte    <J  (0)J  (-t)>,  (1.2) 

0        £    £ 


where  lim  denotes  lim  lim.   Here  the  limit  V  +   »  is  tbe  usual 
thermodynamic  limit,  taken  with  the  density  n  being;  constant;  the 
constants  K  and  the  fluxes  J  ,  which  depend  upon  the  particular  transport 
coefficient  heing  considered,  will  be  defined  in  Chapter  ITT.   The 
brackets  <...>  Henote  an  equilibrium  grand  canonical  ensemble  average. 
(There  are  other  forms  of  these  expressions,  but  all  lead  to  the  same 
final  result.   )   The  fluxes  can  be  separated  into  two  parts:   (1)  a  sum 
over  single  particle  functions  of  the  momentum,  the  'kinetic'  part  .T  ,  and 
(2)  a  sum  over  two  particle  functions  of  the  positions  and  momenta,  the 
'potential'  part  J  . 


k     <j) 

Jk  =  I   a(a) 

a 

a*S  (1.3) 

As  a  result  one  may  separate  three  contributions  to  the  transport 
coefficients  corresponding  to  the  Vinetic-kinetic  term,  the  kinetic- 


potential  cross  terms ,  and  the  potential-potential  terms  in  the  time 
correlation  function: 

?  =  ^  +  2v  +  W  (1-4) 

This  form  is  convenient  for  comparison  with  computer  simulations   ,  and  to 

find  the  first  densitv  correction  only  the  first  two  terms,  £   and  £ 

kk      <j>k 

need  be  considered  since  the  last  term  can  be  shown  to  he  of  higher  order 
in  the  densitv.   In  a  manner  analogous  to  the  cluster  expansion  methods  of 

eauilibrium  statistical  mechanics,  a  density  expansion  of  (1.4)  has  been 

1 2 
derived   ,  allowing  identification  of  both  the  Roltzmann  result  and  first 

density  correction  terms.   ''There  has  been  extensive  discussion  of  the 

resummations  of  this  expansion  which  give  the  non-analytic  portions  of  the 

density  dependence  ;  the  first  two  terms  in  the  series  remain 

unchanged.)   In  Appendix  A,  an  extension  of  this  classical  cluster 

expansion  for  repulsive  potentials  is  given  to  include  quantum  effects  and 

potentials  with  bound  states.   Accordingly  this  work  will  analyze  various 

contributions  to  the  first  densitv  correction  to  transport  coefficients 

for  a  sauare  well  gas  in  an  attempt  to  identify  contributions  due  to  the 

potential  part  of  the  fluxes,  Enskog-like  corrections  to  two  body 

scattering,  three  body  scattering  events,  and  bound  states  contributions — 

in  particular  monomer-dimer  scattering  and  dimer  formation. 

The  results  of  the  following  chapters  are  as  follows.   In  Chapter  IT 

applications  of  the  hard  sphere  fluid  are  reviewed.   The  Enskog  theory   is 

compared  with  molecular  dynamics  and  the  rigorous  results  of  Sengers  et 

al.        The  related  theory  for  time  correlation  functions,  called  the 

1 1 
generalized  Enskoe  theory  ',  is  reviewed  and  calculations  of  the  dynamic 
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structure  factor,  S(k,oj),  are  presenter!.   The  comparison  with  recent 

neutron  scattering  data    and  molecular  dynamics  for  hard  spheres    and 

Lennard-Jones  fluids   is  discussed.   The  modified  Enskog  theory   ,  and  ad 

hoc  adaptation  of  the  Enskog  theory  to  real  fluids  is  also  discussed.   In 

Chapter  III  (and  in  Appendix  A)  the  density  expansion  for  the  hydrodynamic 

transport  coefficients  is  made.   The  quantum  mechanical  form  of  the 

fluxes,  ensemhle  average,  and  cluster  exnansion  is  presented  in  a  compact 

and  general  manner.   The  expressions  giving  the  Roltzmann  transport 

coefficient  and  its  first  order  correction  are  identified 

for  p,  X   and  nl").    The  specialization  to  classical  mechanics  and  the 

sauare  well  potential  is  made  for  t^ese  results  and  is  given  in  Chapter 

IV.   The  Boltzmann  contrihution  is  evaluated  over  a  wide  range  of 

temperatures  and  compared  to  existing  results.   The  contribution  to  the 

first  densitv  correction  bv  £,  is  fjiven  for  in  Chapter  V  for  £  =  n.  A. 

<)>k  ' 

Both  the  scattering  states  and  hound  states  contributions  are  evaluated, 
the  results  agreeing  with  computer  simulation.   Contrasts  are  made  with 
the  behavior  of  real  fluids  and  the  Enskog  theory.   Finally  in  Chanter  VI 
these  results  are  summarized  and  discussed.   The  large  contribution  from 
bound  states  for  £    is  further  analyzed.   Comparison  is  made  between  the 

(pK. 

triple  collision  contribution  to  £.  ,  ,,  based  on  estimates  from  molecular 

kk  ,1 

dynamics,  for  the  square  well  potential  and  Sengers '  hard  sphere 
calculations. 


CHAPTER  II 
HARD  SPHERE  FLUID 

The  hard  sphere  fluid  has  been  a  popular  model  both  because 

theoretical  calculations  are  easier  with  such  a  simple  model  and  because 

for  many  years  the  only  reasonably  successful  model  for  dense  fluids  was 

the  Enskog  model  for  hard  spheres.   The  simple  nature  of  the  interaction 

also  has  led  to  extensive  computer  simulation  over  the  fluid  phase  so  the 

'experimental'  properties  are  well  known,  giving  an  easv  method  of 

checking  theoretical  models.     The  primarv  objection  to  the  hard  sphere 

model  is  an  essentially  trivial  temperature  dependence.   Figure  II . 1 

illustrates  the  difference  between  the  temperature  dependence  of  transport 

coefficients  for  hard  spheres  and  for  the  Lennard-Jones  potential,  in  the 

Poltzmann  Unit.   At  finite  densities  another  qualitative  difference  is 

the  formation  0f  bound  states  (dimerization)  for  potentials  with 

attractive  parts.   Figure  II. 2  shows  the  fraction  of  dimers  for  the 

3 

Lennard-Jones  gas  at  vo~     =  .01.   Thus  it  can  be  argued  that  bound  state 

effects  such  as  nonotner-dimer  scattering  will  have  a  significant  effect  at 
lower  temperatures. 

An  objective  of  this  research  is  to  provide  exact  results  for  the 
first  density  correction  to  the  Boltzmann  results  to  serve  for  guidance 
and  interpretation  of  phenomenological  models  of  the  dense  pas  phase.   As 
mentioned  above,  the  Enskog  theory  describes  dense  gas  transport 
properties  quite  well  over  a  wide  range  of  densities.   However,  the  only 
firm  theoretical  justification  of  the  Enskog  theory  is  given  hv  comparison 
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Figure  II. 1.   RoltzTrann  transport  coefficients  for  the  Lennard- 
Jones  potential  relative  to  hard  spheres.   Horizontal  line  is  hard 
sphere  limit. 
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with  Senders  et  al.  evaluation  of  the  exact  first  density  corrections. 
These  results  and  a  related  calculation  of  the  dynamic  structure  factor 
are  reviewed  brief lv  in  order  to  indicate  the  utility  of  the  Enskog  model 
and  its  range  of  validitv.   The  anplication  of  the  concepts  used  in  the. 
Enskog  theory  to  other  potentials  has  heen  widely  used,  but  prior  to  the 
work  presented  here  there  have  heen  no  exact  results  for  a  critical 
analvsis  of  these  theories.   The  modified  Enskog  theory  is  also  briefly 
described,  and  it  is  shown  to  he  inadequate  for  the  square  well  model. 

o 

A.   Enskog  Theory  of  Hard  Spheres. 

The  Enskog  theory  was  formulated  as  a  phenomenological  modification 
of  the  Foltzmann  transport  equation  using  intuitive  concents  that  are  well 
defined  only  for  a  hard  sphere  fluid.   While  retaining  the  simplicity  of 
two  body  dynamics  Enskog  introduced  modifications  to  account  for  the 
finite  size  of  the  particles,  introducing  the  concents  of  collisional 
transfer  and  exclude^  volume.   The  first  is  essentially  the  transport  of  a 
property  through  the  effect  of  the  interparti cle  potential  rather  than  by 
moving  particles  of  the  fluid  from  one  point  to  another  (Various 
definitions  in  the  literature  of  what  is  meant  by  collisional  transfer  are 
not  exactly  equivalent,  so  care  should  be  exercised  when  comparing 
different  works);  the  second  is  a  modification  of  the  probability  of  a 
collision  due  to  presence  of  other  particles  in  the  fluid.   Unlike  other 
potentials,  the  instantaneous  nature  of  the  hard  sphere  collision  allows 
precise  definition  of  both  the  difference  in  the  positions  of  the 
nartlcles  and  of  the  pair  correlation  function  at  the  time  of  the 
collision.   The  Enskog  equation  gives  the  exact  hinary  collision  dynamics 
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for  hard  spheres  and  correctly  identifies  the  dominant  contribution  from 
higher  order  collisions  to  first  order  in  the  density. 

Comparing  the  theoretical  prediction  to  simulations  for  hard  spheres 
one  notes  the  equation  of  state  is  exact  for  hard  spheres  (see  (2.9)); 

there  is  also  quite  good  agreement  with  the  transport  coefficients  for 

3 

densities  n*  <  . 1  an^  only  w  30%  error  at  n*  =0.3  (n*  =  va~  ,  where  a   is 

the  force  range).   A  summary  of  the  Enskog  predictions  for  the  transport 
coefficients  is  given  below 


nE  =  n0"S'(7(iy  +  *8  ¥   "*  +  '761  *(0)(f-   n*)2) 


H.S.    I        9tt  2ir    2 

xe  =  xo   (77^T  +  1'2  T^  "*  +  -755  gCoXf1  n*)z) 


H.S. 


KR  =  nJ"S'(1.002  g(a)(|1n*)2)  (2.1) 


Here  g(a)  is  the  pair  correlation  function  at  contact  and 

.H.S.    H.S.        H.S. 

A     ,  nn    ,  and  D     are  the  Roltznann  value  for  the  transport 

coefficients : 


H.S.  s         ,m    ,1/2 

0  2        3lT 

16a 


.H.S.    _    ?5   kR      n,      1/2 
64ma 


°  8ma?      B*  (         } 
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Thus  one  sees  that  although  some  density  dependence  is  included,  the 
temperature  dependence  is  the  sane  as  in  the  Roltzmann  case.   A  comparison 
of  the  Enskog  predictions  for  transport  data  from  molecular  dvnanics1 '  is 
shown  in  Tahle  II. 1. 

Direct  comparison  of  the  Enskog  theory  to  experiment  does  not  show 
such  good  agreement;  however,  hy  using  the  low  density  transport  data  to 
define  a  temperature  dependent  hard  sphere  diameter,  one  gets  reasonably 
good  agreement  with  experiment"   (st  107  error  at  the  critical  density). 

B.   Sengers  et  al.  Calculation  of  the  First  Density  Correction  for  Hard 
Spheres 

As  mentioned  above,  the  simple  nature  of  hard  sphere  dynamics  makes 

theoretical  calculations  easier.   In  particular  only  sequences  of  binary 

collisions  occur  so  that  the  binary  collision  expansion    of  the  kinetic 

equation  for  hard  spheres  is  particularly  useful.   It  can  be  shown  that 

for  three  hard  spheres  interacting  in  three  dimensions  there  is  only  a 

finite  number  of  binary  collisions;  so  the  expansion  terminates  after  a 

few  terms.   Using  this,  Sengers'  group  was  able  to  evaluate  the  first 

density  correction  to  the  hard  sphere  transport  coefficients.    The 

results  are  in  good  agreement  with  the  Enskog  values  (r  107),  and  they 

were  able  to  identify  which  types  of  collisions  were  included  in  the 

Enskog  theorv  and  which  were  neglected.   A  parallel  of  this  analysis  with 

the  square  well  model  is  not  possible  exactly;  however,  some  of  the 

collision  sequences  for  scattering  by  scmare  wells  can  be  related  to 

similar  sequences  for  bard  spheres  and  may  thus  be  considered  as  having 

the  same  relative  imnortance  for  determining  the  square  well  transport 

coefficients. 
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TABLE  IT.]  Comparison  of  Enskog  Theory  to  Hard 
Cohere  Molecular  Dynamics 


no  D/nRa       D/nFb         n/n^,        A/X^       k/<e 

.0071  1.00  1.02  1.01        .98  .9 

.0354  1.01  1.04         1.01       .99         1.0 

.0707  1.03  —  .99 

.1414  1.00  i#lf;  .oa         .07  #q 

.1768  1.13 

.2357  1.16  1.34         1.02      1.00  .98 

.2828  1.17 

.3536  1.06 

.3928  .95 

.4419  .70 

.4714 


Note:   Ratio  of  hard    sphere  transport  coefficients  to  the  values  predicted 
by  the  Enskog  theory.   Column  headed  hv  D/D    is  for  108  particle 
simulation;  column  headed  by  D/D    is  an  extrapolation  to  infinite  system, 
using  mainly  a  500  particle  simulation.   Ouoted  uncertainities  are 
*  ±  .02  for  D/DE,  ci  ±  .04  for  t\/t\        a  ±  .02  for  A/A   and 
si  ±  .02  for  k/k.  E 

Source:   B.  J.  Alder,  D.  M.  Cass,  and  T.  E.  Wainwright . *° 
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Here  the  density  expansion  is  ^  =  F   +  (no'  )g  .   The  first  terms  in 
X^  and  n   are  due  to  the  potential  parts  of  the  fluxes  and  are  given 
exactlv  hy  the  Enskog  theory. 


C .   Generalized  Enskog  Equation 

A  useful  property  of  hard  spheres  is  the  Instantaneous  nature  of  the 
collision.   This  implies  that  one  can  calculate  in  detail  the  t  ■*■  0+ 
limit  in  the  derivation  of  a  kinetic  equation  and  still  retain  collisional 
effects  associated  with  transfer  of  momentum.   This  procedure  leads  to 
what  is  called  the  generalized  Enskog  equation.   (For  continuous 
potentials,  it  leads  to  a  collisionless  Vlasov-like  equation.)   This  short 
time  limit  does  not  require  a  density  expansion,  and  therefore  the 
applicability  of  the  kinetic  equation  is  not  limited  to  dilute  erases.   As 
an  example  the  dynamic  structure  factor,  S(k,u),  determined  from  the 
p-eneralized  Enskog  theorv  (GET)  shall  be  reviewed.2^ 
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S(k,oj),  which  is  related  to  the  differentia]  cross  section,  is  given 
by  the  Fourier  transform  in  space  and  tine  of  the  density-density 
correlation  function, 

S(k,o>)  =1/"   dteiU,t  /d?e-^,?<nro,n);c?,-t)>, 

—CO 

(2. A) 

with 

*  +  N 

n(r,-t)  =   I   6(r"-q\(-t)). 

i=l 

Following  the  usual  methods  of  kinetic  theory,  one  notes  that  the  solution 
of  (2.4)  requires  a  forma]  closure  of  the  RBGKY  hierarchy  hv  expressing 
the  two-body  reduced  distribution  function  as  a  functional  of  the  one-body 
reduced  distribution  function  (see  Appendix  A).   For  the  type  of 
correlation  function  considered  here,  suns  of  single  particle  functions, 
this  is  relatively  simple  because  there  exists  a  set  of  operators 
U(l...s;t)  which  map  the  single  particle  function  onto  the  reduced 
distribution  function,  (i.e.,  i|;(l...s;t)  =  U(  1.  .  .s;  t)6(r-q  )  for  this 
case.)   Thus  one  can  easily  close  the  hierarchy  to  get  the  kinetic 
equation  for  i|»(l;t): 

(j£   +  ~  '   ^-)lp(1)(l;t)  =  -n/dx?TJT(12)H(12;t)Trl(l;t)  *(1)(l;t). 

(2.5) 

Here  Lj  is  the  interacting  part  of  the  Liouville  operator. 

rT  ,,,,>.      3  d>C  12)    3 

I  L,  yi/-)    = •  — - — ■ tor  continuous  potentials,  and 

1  8rl2,i     ^12, i 

Tjj(12)  =  b^   tnp  binary  collision  operator  for  hard  spheres.   ]   The 

usual  density  expansion  is  obtained  from  a  cluster  expansion  of 
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tJ(12;t)Tl   (l;t);  however  by  taking  the  t  +   0   limit  of  U(12;t)U-1(l  ;t) , 
one  gets  the  generalized  Enskog  equation.   As  in  the  Enskog  theory  this 
procedure  retains  static  correlation  between  the  particles  and  binary 
collision  dynamics,  and  thus  one  can  consider  the  GET  as  giving  a  valid 
approximation  to  the  time  correlation  functions  over  the  same  range  of 
densities  as  the  Enskog  theory.   "For  S(k,a>)  the  Fourier  transformed  form 
of  (2.4)  leads  to  the  following  expression: 

S(k,u))/S(k)  =  ?  Re  (1,  R(k,u)  L).  (2. ft) 


The  scalar  product  is  defined  with  a  Maxwell-Boltzmann  weighting  functi 
f0(v): 


(a,b)  =  /  dv  f  (v)  a*b,  (2.7: 


and  the  operator  R(k,uj)  by 


R(k,o>)  =  [-iu)  +  ik«v  -  R(k)]    ,  (2.8) 


with  B(k)  being  the  generalized  Enskog  operator  defined  by  equation  (2.5) 
of  reference  (24).   (The  transport  coefficients  obtained  from  (2.8)  are 
the  same  as  those  from  the  Enskog  theory.)   One  notes  that  B(k),  derived 
from  the  right  hand  side  of  ''2.5),  does  not  depend  on  to  since  the  t  ■*■  0 
Unit  was  taken,  but  it  still  generates  dynamics  due  to  the  discontinuous 
potential. 

Due  to  the  complicated  form  of  P(k),  one  usually  proceeds  by  choosing 
a  kinetic  model  for  R(k,u).   This  involves  finding  the  matrix  elements  of 


Figure  II. 3.   Dynamic  structure  factor  S(k,u))  vs.  k  for  to  =   0,1,2 
meV.   Kinetic  theory  (-s— ) ;  neutron  scattering  data  for  Krypton 
(ooo)  at  T  =  297°K,   no  =  .47. 
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Figure  II. 4.   Full  width  at  half  maximum  (FWHM)  vs.  k.   Kinetic 
theory  (— — );  Krypton  data  (ooo);  Lennard-Jones  molecular 
dynamics  (•••);  hard  sphere  molecular  dynamics  (AAA). 
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R  with  respect  to  a  finite  set  of  orthogonal  functions,  |u  },  which 

guarantee  that  some  set  of  desired  properties  of  the  correlation  function 

is  given,  while  representing  the  hehavior  outside  the  space  spanned  hy 

ju  }  in  terms  of  a  single  parameter.   By  choosing  the  \u    }    so  the 
L  a'  l  or 

bydrodynami c  limit  of  S(k,u)^  is  given  correctly,  very  eood  agreement  is 
found  with  the  molecular  dynamic  results  of  Alder  et  al.  ""  for  hard 
snheres,  and  fair  agreement  is  found  with  the  neutron  scattering  data  of 
Eglestaff  et  al.     These  results  are  shown  in  figures  II. 3  and  II. 4,  and 
one  sees  that  there  are  noticeahle  discrepancies  hetween  the  Krynton  data 
and  the  GET  at  snail  k  and  to  values.   Included  in  figure  II. A  are  the 
molecular  dvnamic  results  of  Ullo  and  Yip  for  the  Lennard-Jones 
potential.     The  full  Lennard-Jones  model  is  closer  to  the  experimental 
data  over  a  larger  range  of  k  values  although  some  differences  are  seen  at 
higher  density  in  the  small  k  range  (k  <  1  A   ).   The  truncated  Lennard- 
Jones  potential  (no  attractive  forces)  behaves  more  like  the  hard  sphere 
kinetic  mode]  at  small  k;  therefore,  the  differences  hetween  the  OFT  and 
experiment  are  due  to  the  potential  model  not  containing  attractive 
forces. 


n.   Modified  Fnskog  Theory. 

From  the  results  of  the  Fnskog  theory  one  can  make  various 
phenomenological  modifications  for  real  fluids.   por  example  if  one 
assumes  that  the  processes  the  Fnskog  theory  includes  are  the  important 
ones  for  real  fluids  as  well,  then  one  can  claim  the  form  of  the  transport 
coefficients  is  the  same,  with  the  low  densitv  result  and  g(oO  to  be  found 
from  experimental  data.   The  pair  correlation  function  at  contact  cannot 
be  unambiguously  defined  for  other  potentials,  however,  an  ad  hoc  method 


28 


of  defining  g  through  the  equation  of  state  has  heen  extensively  tested. 
Ope  knows  that  for  bard  spheres  the  equation  of  state  is  related  to  g(a) 
by 


£  [(*g£) -»]--?  «c>-.  »•» 


?y   using  experimental  data  on  the  left  hand  side,  Hanley  et  al.    have 
made  comparisons  of  the  predictions  of  the  modified  Enskog  theory  (MET) 
with  experiment.   Their  results  show  the  MET  agrees  quite  well  with 
experiment  at  low  densities  over  a  wide  range  of  temperature  and  has  an 
error  of  ^10%   at  densities  one  or  two  times  the  critical  density.   They 
also  find  that  the  first  density  correction  (5,,  if  5  =  Kn   +  5, n  +  . . . ) 
predicted  by  the  MET  agrees  with  experiment  if  experimental  virial 
coefficients  and  low  density  transport  data  are   used.   On  the  other  hand 
if  the  virial  coefficients  and  Roltzmann  transport  coefficient  for  a 
Lennard-Jones  potential  are  used,  then  poor  agreement  is  found.   This  is 
attributed  to  the  failure  of  the  MET  for  the  Lennard-Jones  fluid  to 
predict  the  triple  collision  contribution  for  real  fluids  to  E, 
correctlv.   Presumably  then,  computer  simulations  for  a  model  fluid  would 
he  well  given  by  the  MET  for  that  model.   However,  Micbels  and 
Trappeniers   have  connared  their  molecular  dynamics  simulation  for  the 
square  well  fluid  to  tbe  MET,  and  they  conclude  that  it  does  not  give  even 
the  correct  sign  of  the  first  density  correction.   It  is  not  known  if  the 
same  result  holds  for  continuous  potentials  such  as  the  Lennard-Jones 
model . 

In  figure  II. 5  the  MET  value  for  the  first  density  correction  to  nD 
for  the  square  well  is  compared  to  molecular  dvnamics  and  the  predictions 


Figure  II. 5.   Comparison  of  modified  Enskog  theory  (-  -  -), 
repeated  ring  theory  (       ,  with  hound  states;  -  o  -  o  - 
without  bound  states),  and  molecular  dynamics  (A)  for  the  first 
density  correction  to  self  diffusion  for  the  square  well. 
Horizontal  dashed  line  is  the  hard  sphere  limit. 
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2  1 
of  a  'repeated-ring'  model."    It  is  noted  that  the  MET  result  is  similar 

to  that  of  Hoffman  and  Curtiss  ",  while  the  nrimarv  difference  hetween  the 

Hoff man-Curt iss  model  and  the  repeated  ring  calculation  is  that  the 

repeated  ring  pronerly  includes  the  effects  of  hounds  states,   Thus  it  nay 

he  seen  that  boun^  states  or  contributions  due  the  attractive  part  of  the 

well  are  important  in  dense  ^ases. 

From  the  examples  ahove  one  can  see  that  the  hard  sphere  model  can 

give  qualitatively  correct  hehavior  for  dynamic  properties  of  gases,  and 

the  agreement  hetween  hard  sphere  molecular  dynamics  and  theory  is 

generally  verv  good.   However,  the  S(k,u)  calculation  and  the  failure  of 

the  square  well  MET  results  show  that  attractive  forces  and  bound  states 

effects  play  an  important  role  in  the  calculation  of  dynamic  properties. 

An  alternative  dense  gas  theory  must  then  take  into  account  effects  which 

are  not  included  in  or  properly  calculated  by  models  hased  on  hard 

spheres. 


CHAPTFR  III 
FORMAL  DENSITY  EXPANSION  OF  THE  TRANSPORT  COFFFICIENT 

The  densitv  expansion  of  time  correlation  function  for  the  transport 

coefficients  has  been  discussed  previously  for  classical  mechanical 

systems  by  several  authors.  '      The  generalization  given  here  for  quantum 

mechanical  systems  follows  in  an  obvious  way  from  the  classical  results; 
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McLennan   and  Boercker  and  Dufty  "   have  considered  the  special  case  of 

correlations  between  sums  of  single  particle  functions  for  quantum  systems 

previously.   The  transport  coefficients  considered  here  are  of  the  general 

form 


A  =  Pp  lim  -P-  j     dte    <J.  •  J,(-t)>,  (3.1) 

V  J  A    A     ' 


where  K.    is  a  constant  depending  upon  the  particular  transport  coefficient 

considered  and  the  brackets  <...>  denote  a  grand  canonical  ensemble 

average.   The  transport  coefficients  A  for  a  simple  fluid  are  the  shear 

viscosity  (p),  the  thermal  conductivity  (X),  the  hulk  viscosity  (k),  and 

self  diffusion  (nD,  where  n  is  the  number  density).   In  the  definitions  of 

J,  given  below  and  in  the  correlation  function  J,  •  J.(-t),  it  is  assumed 
A  A    A    ' 

that  repeated  Latin  indices  are  summed  over  (i,i,  etc.  =  x,y,z): 

Self-diffusion:   K  =  ^ 

3m  " 


(3.2a) 


Jn  =  /N  Ph   > 
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Shear  viscosity:   K   =  — - 

n   10  (3.2b) 


Jn^  lT±1>Tr» 


V 
Thermal  conductivity:   K,  = 


X  3 


J   =  S   -  -  P 
A    i    pi' 


(3.2c) 


Bulk  viscosity:   K   =  Tr- 
IC     9 


(3. 2d) 


J  =  T .  .  -  f   ' 

K      11      11 


In  (3. 2d)  f        is  the  time  average  of  T,  ,  .   McLennan  has  shown  T   is  given 
kk  k.k  kk 


by 


?kk  =  3^V  +  (!&,<«  ~   <H»  +  (&(N  -  <N»).  (3.3) 


kk    *•      ^3u;n  ^3n 


In  (3.2)  and  (3.3),  in  is  the  particle  mass,  N  is  the  numher  of 
particles,  3  =  — - -  with  k„  being  P.oltzmann 's  constant  and  T  the 
temperature,  p  is  the  mass  density,  h  is  the  enthalpy  density  (h  =  u  +  p, 
with  u  the  energy  density  and  p  the  pressure),  and  F  is  the  Hamiltonian 
operator.   The  braces  {••},„   denote  the  traceless  part  of  T..: 

6  . 

1   ;Tr    ii     3   kk'  •J*^ 


with  6.  .  beine  the  Kronecker  de] ta  function.   The  stress  tensor  T.  ., 
ii  l.l ' 

energy  flux  S.  and  total  momentum  P.  are  defined  by 
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1   v   3*(o8) 


t    =  -  y  -  -  y 

a  a*g   ag,i 


1        2       1 

s-  =  — o     I   D  P  ■  +  7—  I  U(ag)(p  .  +  p„.  ) 
i   0  2  L     a  ax   4n  f;0  lY   ^   ai    gi 


-  (p  .  +  PQ.)  -^r — —  r  .  .  (3.5) 

ai    gi   9r      ag,i J  w   ; 

ag,  i 


p.  =  y  p 

1    L      c 


Here  a,g  are  particle  labels,  p   is  the  momentum  operator,  r   =  r   -  r 

a  ag    a    g 

is  the  relative  Position  operator,  and  <j>(a,g)  in  the  pair  potential.   Th< 
operation  indicated  by  [a, hi,  denotes  the  anti-commutator  ah  +  ha.   Tt  is 

-r 

seen  that  the  fluxes  in  general  have  two  parts:   ( 1  )  a  'kinetic'  part 

consisting  of  sums  of  single  particle  function  of  the  momentum,  J  ;  (2)  i 

'potential'  part  consisting  of  sums  of  two  particle  functions  of  the 

momenta  and  relative  position,  J  ;  i.e., 

<f> 

J   =  I   a(o) 


J  ={  I     b(ag) 
*     a*g 


There  exist  several  other  definitions  for  the  quantum  formulation  of  the 
time  correlation  functions  for  the  transport  coefficients;  however,  they 
are  all  equivalent  for  the  transport  coefficients  considered  here. 
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From  the  definition  of  J  one  can  write  the  transport  coefficient  as  a 

sum  of  three  terms, 

A  =  A.,  +  2A..  +  A.  .  ,  (3.6) 

kk     <j>k    <(>9 

caller!  the  kinetic,  cross,  and  potential  terms,  respectively.   This  form 

is  convenient  for  comparison  with  molecular  dynamics  since  each  term  in 

(3.6)  can  he  calculated  senarately.   It  will  he  shown  below  that  A.,  is  of 

kk 

0 
lowest  order  in  the  density  (leading  tern  •>  n  )  and  gives  the  usual 

Boltzmann  result.   The  first  density  correction  to  A  (order  a  n  )  then 

consists  of  the  second  term  in  the  expansion  of  A    plus  the  lowest  order 

term  in  A,,  ,  while  A,  ,  is  formally  of  hi<rher  order  in  n.   The  contribution 
<j>k         $*)> 

from  A    depends  upon  two  body  dynamics  only  and  will  be  evaluated  in 
<pk 

Chapter  IV.   The  contribution  from  A,,  involves  three  body  dynamics. 

kk 

The  details  of  the  density  expansion  of  each  term  in  (3.6)  are  given 
in  Appendix  A.   Compact  and  general  expressions  can  be  found  if  the 
statistics  of  the  system  are  included  in  the  density  operator,  thus  the 
reduced  distribution  function  and  cluster  functions  can  he  defined  in  a 

97 

manner  closely  related  to  the  classical  calculations."    It  is  to  be 

note'1  that  the  cluster  expansion  of  the  densitv  operator  then  implies  that 

3  3 

both  no   (where  a  characterizes  the  force  range)  and  nx"  (where 

8   1/2 
x  =  h(— — )    is  the  thermal  wavelength)  are  small  parameters.   This  grand 
2irm 

canonical  average  will  therefore  have  the  form 


N 
N>0   '"   1..N 


-   3N   -^ 
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N 
3  =  Z  ~  Tr  p   .  (3.7) 

N>0  '  '  1..N 
Here  the  trace  is  taken  over  a  complete  set  of  states  in  the  N-particle 

space,  while  the  operator  P   projects  onto  the  properly  symmetrized 

3     ftu 
suhspace  for  identical  particles;  the  activity  is  defined  by  A  z  =  e   , 

where  u  is  the  chemical  potential.   In  a  manner  analogous  to  the  classical 

case,  reduced  distri hntion  functions  corresponding  to  the  kinetic  and 

potential  parts  of  the  fluxes  can  he  found,  and  each  of  the  terms  in  (3.6) 

can  be  represented  as  a  trace  over  the  coordinates  of  one  and  two 

particles  involving  the  one-  and  two— bodv  reduced  distribution  functions: 


A   =  Re  lim  ^  Tr  a(l)  *£  (l;e) 

2 

A   =  Re  lim  ^  Tr   b(] .2)tJ2)(12;  e) 

<()k  IV    l    0  k 

2  .    . 

h.±   =  Ee  lim  ^-  Tr   b(12)\K  ;(12;e)  (3.8) 

99  lv    x    2       * 

( s )      ( s ) 
Here  ty  and  ip    are  the  kinetic  and  potential  parts  of  the  s-particle 

k  a) 

(s) 

reduced  distribution  function,  lj;    ,  defined  bv 


+(8)(l..S;e)  =  /  dte-£t$Cs)(l..s;t) 

n 

(    ^  N 

M>s      '     S+1...N 

(1)    (2)       (2) 

The  density  expansions  of  ik   ,  ik   ,  and  ill    therefore  give  the  density 

k     k         cp 

expansion  for  A.   The  results  of  Appendix  A  imply  for  each  term  in  (3. ft) 
that 
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\k  =  \k,n  +  n  V.i 


A.,  =  0  +  P  A     +  ... 


A„  =  °  +  n  +  "\*   7  +  •••  (3.10) 


(?> 
Thus  t^e  cluster  expansion  for  ty   '       is  not  needed  for  the  first  density 

correction.   Also  as  in  the  classical  case,  ill,    is  determined  from  the 

k 

first  equation  of  the  BRGKY  hierarchy: 


(e  +  h(l))^])(l)  +  n  Tr  L];(lv)T(i2;e|^1)(.)>  =  ^1)(l;t-0).   (3.11) 


Here  T  is  the  expression  for  ik    as  a  functional  of  U»,    as  piven  hy  the 

k  rk 

cluster  expansion, 


/2)(12;~ ^ 


%    '(12;e)  =  ^(12;e|^;J' (•)).  (3.12) 


The  Liouville  operator  in  (3.11)  is  defined  by 


L  =  I   L(a)  +  I     L  (a&)  =  £  [H,  ]_  ,  (3.13) 

a        a<g 

where  [A,B]_  is  the  commutator  of  A  and  B,  AB  -  BA.   Since  "F  is  a  linear 

operator  on  ik   ,  (3.11)  may  be  formally  solved  for  ip    : 


(1 )  "   -1  «ni 

<JV   (l;e)  =  (e  +  L(l)  +  m')   iu;(l;t=0) 
k  k 
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M'f(l)  =  n  Tr  L  ( 12)*(12; e I f ( • ) ). 


(3.14) 
Pirect  suhstitutiop  of  the  density  expansion  of  ¥  from  Appendix  A  gives 


the  density  expansion  for  M*.   Noting  that  only  the  e  ■»•  0  limit  of  (3.14) 
is  required  in  (3.8)  or  (3.10)  and  that  i^    is  independent  of  r,  , 
equation  (3.14)  can  he  rewritten  as 


g(l)  =  M  ^(1), 


(3.15a) 


with 


t(l)  =  lim  n  ](l)  n^(1)(l  ;e), 
e+o         k 


h(l)  =  n  Vl)  f(1)(1;t=0), 


_] 
and  M  =  lim  D   (l)M'P(l). 

e+0 


(3.15b) 


3  — fiHf 1 ^  A 

Here  D(l)  =  X  e       .  The  density  expansion  of  M  (from  f)  and  h( 1 )  (from 

a  cluster  expansion  of  (3.°))  identifies  the  density  expansion  for  g(l). 

Letting 


g(l)  =  g(0)(l)  +  n  g(1)(l)  + 


h(l)  =  >/0)(n  +  n  ha)(l)  + 


Ml  +  n  M2 


(3.16) 


then  the  lowest  order  terms  in  (3.10)  will  he  given  hy 
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A..     _  =   Re    lim  ^  Tr   D(l)    a(n)(l)    g(0)(l) 

kk,0  „        V    , 


A,  ,     ,    =   Re    lim  -^  Tr   D(l)    (a(1)(l)    g(0)(l)    +  a(0)(l)    e(1)(l)) 
Ick  ,1  V 

'  V+oo        1 


A..  .  =  ^    Jtm-^Tr   b(0)(12)  lAl2;\v{l)gC0)(l)), 
♦k'J      V—  2V  1,2  * 


(3.17) 


(0)   (i)      en) 

Here  a    ,  a    ,  and  b    are  the  lowest  order  terns  in  the  density 
expansion  of  the  fluxes  (this  density  dependence  is  only  through  thermo- 
dynamic functions),  and  4*   is  the  first  term  in  the  density  expansion 
of  V  [specifically,  f  (12|f)  -  lim  r>(  12)eR(  12  ;  e)(l+P   )P-1(  1  )f  (  1 )  ]  , 
equation  (A. 27).   The  solutions  to 


VC0)(l)  =h(0)(i) 


and 


M]g(1)(l)  =-N2,(0)(l)  +hCn(l) 


(3. 18) 


Rive  g-   (I)  and  g    (1).   Defining  the  scalar  product  (a,h)  allows  a  more 
compact  expression  for  the  kinetic  parts  of  (3.16): 


(a,h)  =   lim  j  Tr  D(l)a  (l)b(l) 

V-*-°°   1 


(3.19) 


where  the  dagger,  t,  denotes  the  hermitian  adioint.   Thus 


A,  ,  _  =  Re  (a(0)(l),g(0)(l)) 
kk,0 


JL,  ,  =Re  [(a(n)(l),g(1)(l))  +  (a(1)(l  )  ,g(0)  (1)  ) 

KK  ,  1         L 
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AAI     ,    =   Re  ^77  Tr      h(0)(]2)    ¥(12|D(1)R(0)(1))    . 
<I>V>1  2V        0 


(3.20) 


Using  the  definitions  (3.2)  of  each  transport  coefficient  in  (3.20)  gives 

the  specific  equations  which  must  he  evaluated. 

For  the  shear  viscosity,  n  =  n,  .  n   +  n(n,  ,  ,  +  2n,  ,)+...: 

kk  ,  0  kk  ,1  <p  k  ,  1 


'kk,0      10   v    ii'    ii    ' 


%v    ,    =  7n    CD..,10)) 

kk,  1         10        n      ij 


^     i    =  T77F     Tr    [<M12>ri9    ,P,,    •)    *    *,  (12|D(1)    fi(?})] 
tp^jl         10V    .    .    L      I  12,1    12,i  1  ]i       J 

i ,  z 


(3.21) 


Here 


^PliPl^Tr    ' 


and 


8(0)  -1 

P. .        =M,      D. .    , 
il  1        ij 


^--m:1  kk1  d. 


1        2    1        ij 


(3.22) 


For  the  thermal  conductivity,  X  =    X,  ,  .  +  n(X.,  ,  +  2X,,  ,)  +  ...  : 

kk  ,0      kk  ,  1      <j)k ,  1 


kk ,  0    3 


'S<0)(1),^0)(1) 


K, 


V 


kk ,  1     3 


(s!0)(i),§a)(D) 
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V,i   =  ~  f2  ^  t*(12)'(pli  +p21)]+ 


+  ^[[(P!j  +P2i)'  Li(12)Pl2,il>  ^^U 


(3.23) 

f.(l,2|n(] )§f0)(l))}. 
1  i      J 


Here 


5(0).n  _  1       "1    ,(0) 
S.   (1  )  =  —  p   .  (- h    ) 

l         n   1,1   2m 


Mn)(n=M^s;n)(i) 


3(0)(1)  =  M-^M"3  S<0)(1).  (3.24) 


Note  that  tv  '  is  the  first  term  in  the  density  expansion  for  the  enthalpy 
per  particle.  One  might  expect  there  would  he  a  tern  involving  the  second 
order  term  in  the  enthalpy,  h  ,  however  it  vanishes  due  to  orthogonality 
conditions. 

Finallv,  for  the  self-diffusion  coefficient,  (nD)  =  (nT1)n 
+  n(nn),  +  .. . : 


(nD)Q  =  -^   (Pl.,^0)(l)) 


3m 

(ND)j  =-^  (p11,?J1)(l))  (3.2S) 

3n  ' 


with 
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JJ(O) 


CD  =  m  !P,. 

si   ll 


P    (1)  =  M  !m  M  Jp    .  (3.26) 

si  s2  slh  w    ' 

Here  Mg  is  the  collision  operator  for  self  correlations  defined  in 
Appendix  A. 

The  classical  Unit  of  these  equations  for  r|,  A  and  nD  will  be 
considered  in  Charter  IV.   The  classical  limit  agrees  with  the  results  of 
Ernst  et  al.   ,  and  the  quantum  operator  M~  agrees  with  that  of 


McLennan 
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CHAPTFR  IV 
CLASSICAL  SOUARE  WELL  CAS 

The  square  well  potential  defined  by  (1.1)  simplifies  the  calculation 

of  the  transport  coefficients  since  (!)  the  time  evolution  of  the  system 

can  be  written  in  a  closed  form,  and  (2)  the  discontinuous  nature  of  the 

potential  reduces  the  initial  conditions  to  he  considered  to  surfaces  at 

a  and  Ro.   From  Appendix  B  and  equations  (A  .26)  the  classical  expressions 

for  the  collision  operators,  M   and  M   ,  may  be  written  in  the  following 

forms : 


M  g(l)   =  -   lim   f01(l)/dp2dr   8(12)eR(12)    f    (12)(l+ffP      )    g(l). 
e+0 


M2g(l)    =  fn1(l)/dp2dp3dr'21dr*31    6(12)[eR(123)W(123)    -  eR(12)W(12) 


x   (eR(13)V(13)   +  eR(23)W(23)   -   l)(l+oP     +aP      )f   (2)f    (3)r(1).        (4.1) 


The  classical  W-functions  are  W(l...s)  =  exp(~3  .E.  <(>(ii))i  ^nr!  the 

i<l 

..       3<t>(r)  3    r  .  .,..,-,, 

0-onerator  is  given  by  ~— '■ - —  for  continuous  potentials;  it  will  be 

3r.    3p. 
l     i 

assumed  that  the  potential  is  continuous  at  this  point,  and  later  the 
limit  to  a  square  well  potential  will  be  used.   Then  the  transport 
coefficients  may  be  written  as  the  following:   (1)  For  the  shear 
viscosity: 
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\k,o=To/dPi  fo(1)Di.i(1)"x.i(1) 

\k,l    ^/d?l    V1)fii.i(1)   M2Bij(1) 


-»--»-•*- 


^i     i    =  ~  7n   lirn  /dp,dp_dr   (r,e(12)p.)eR(12;e) 
d>K,.l  40        „    '       12  i 


x   f    (12)(l+P1,)fl.,(l).  (4.2) 

0  1 2      1 1 


Here  f^Cl)  and  f,-,(12)  are  the  low  density  one-  and  two-body 
distribution  functions: 


f0(12)  =  V,)f0(2)  e"e*(12) 


v 

fnCl)  -  (f  )"3/2  e"B  ?m   ,  (4.3) 


and  r  =  r   -  r   and  p  =  — -  (p   -  p  )  are  the  relative  position  and 

mo™entu"i.   D..(l)  =  lp,,p,.}m   is  the  classical  kinetic  part  of  the 
1.1       l  li  1.1  '  Tr 

momentum  flux;  D   (1)  is  the  solution  to 


M  f)       =  D.   ,  (4.4) 

1  i.i    ii 


where  M  is  seen  to  reduce  to  the  classical  Roltzmann  operator.   The 
potential  part  of  the  momentum  flux  has  been  written  as 


:.  I*-  £  r.  9(12)p.  .  (4.5) 

l  dr,     l        i 
.1 
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This  allows  the  cross  term  to  be  written  in  the  su^pestive  form 


^  i  =  "777  lim   /dp1dp,dr(r.p.e9(12)R(12)fn(12)(l+P..)B..(l),    (4.6) 
<pk ,  1    40    „  '   1   2     ii  (i         1 2   1 1 


where  the  combination  of  operators  lim  e9(12)R(12)  can  he  evaluated  in  the 
same  manner  as  in  the  Boltzmann  case  (see  Chanter  V).   The  solutions  to 
(4.1)  depend  upon  the  evaluation  of  P . . .   The  solution  to  equation  (4.4) 

for  D..  is  usually  £iven  in  terms  of  successive  Sonine  polynomial 

12 

approximations   ;  here  the  first  Sonine  polynomial  approximation  which  is 

known  to  accurate  to  within  a  few  percent  will  he  used;  i.e., 


B,  .  -  a  nD.  .  (4.7) 

ij    n,0  n 


where  this  in  turn  implies  from  (4.2)  and  (4.4) 


\k,o  =  g"lan,o 


where    a      „   is    siven   hy    (^.1Q)    helow.      For    the    other    transport    coefficients 
n  ,0 

similar   results    follow: 

(2)    For    the    thermal    conductivity: 


K^/pV  s1(i).^1(i) 


kPB      r.+   .    ,+   .M 


kk ,  0  3 

2 

xkk,i  =-T^/^1V?i)§i(1),M2§ia) 

? 
k   3" 

X..     .    =   lim  -5-5-  Jdp.dp   d?  «j,(r)P,fn(12).eR(12)(l+P      )5 ,(1) 

♦k»!        e+0    12m'  :  x   n  12      x 
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2 
kg 

+  lim  -2—  /dp  dp  dr  r.(P»p)8(12)»eR(12)f  (12)(1+P  _)§  (1).     (4.8) 

e+0  12m  X 


Here  S.  is  the  low  density  form  of  the  kinetic  part  of  the  heat  flux 

2 

V1)  =  ^r-T3)pii  ■  (4-0) 


md   S.(l)  is  given  by 


§  (1)  =  >f~1S.  (1).  (4.10) 


In  the  first  Son^ne  nolynorial  approximation,  S   is  given  by 


3,(1)  -a^S^l) 


where  the  first  of  equation  (4.6)  and  eauation  (4.10)  then  give 


5k 
\k,0  =  ^mj  aX,0'  r4']1) 


when  a    is  given  by  (4.19).   ^gain  the  potential  part  of  the  flux  that 
depends  upon  the  force,  usually  written  in  the  form  (P»  ( 3<j>(r )  /3r)  has  been 
eliminated  in  favor  of  9(12)(P'p)r,  . 
(3)  For  self-diffusion: 


(nD)0  -— J  /dPl  V1)pli#?li 


(nD)1  =fj!   ^    f0(D?11-M2s?u   .  (4.12) 

Jm 
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Here  P   is  given  by  the  inverse  of  the  Lorentz-Boltzmann  operator, 

M    ,  acting  on  p,  ,  and  M   in  the  triple  collision  ooerator  for  self 
is  J-       z  s 

correlations , 

?li  =  ^s  Pli  '  ^'13> 

M   and  M    are  given  by  (4.3)  and  (4.4)  with  a  =   0.   In  the  first  Sonine 
polynomial  approximation, 


11 

with 


(nD)   =  (n3)  *a  ,  (4.14) 


where  a^  is  define  in  (4.19)  below.   Equation  (4.2),  (4.8)  and  (4.12) 
for  the  Boltzmann  term  and  cross  term  depend  upon  the  evaluation  of  the 
operators  eR(12)  and  9(12)eP(12)  acting  on  single  particle  functions  of 
the  momenta.   For  scattering  states,  it  is  easy  to  see  that  the 

lim  eR(12)fj(l)  =  lim  g(l;-t)  if  g(-t)  ■»■  a  constant  as  t  +■  ».   For  bound 
states  this  identity  no  longer  holds,  an^  more  care  must  be  taken  in  the 
evaluation  of  the  limit  e  +   0.   In  addition,  the  0(ii)  operator  for  the 
square  well  is  undefined  and  must  be  evaluated  in  the  context  of  a  limit 
of  a  continuous  potential.   The  presence  of  the  9  operator  restricts  the 
integral  to  those  portions  of  phase  space  where  the  force  is  non-zero,  and 
since  the  9  operator  always  occurs  in  conjunction  with  the  resolvent  R, 
the  following  identity  can  be  used: 


)(ii)R(ii)  =  1  +  R  '(ij)R(i.l).  (4.15) 
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The  use  of  (4.15)  and  Gauss'  theorem  reduces  the  integration  over  the 

relative  coordinate  to  the  surface  hounding  the  space  -r-*- — —  #  0,  plus  some 

or . 

integrals  which  vanish  or  e  -*-  0,    i.e.,  integrals  of  the  forin 


lim  /dp  dp  /dr  b(12)e9(12)R(l2)g(12) 

=  lim  /dp  dp  /  dr  h(12)  e  (1  -  R_1(12)R(12))g(12) 

e+0         V 

s 

=  lim  -  /dp  dp\  /  dr  h(12)(-  p.  -|-  )  eR(12)g(12) 
n  1   2  ',  *"m   l  dr. 


e-i-n 


lim  -  /dpdp0  /  dr  lr|2h(12)  -(p»r)eR(12)g(12) ,       (4.16) 


e+o        s 


3  (J» 


where  V   is  the  volume  where  _    ^  D  and  S  is  the  surface  bounding  that 
s  dr. 

l 

volume.   Thus  one  can  see  that  the  square  well  integrals  for  the  Boltzmann 
and  cross  terms  can  be  expressed  in  terms  of  integrals  over  the  surfaces 
|rl  =   a    ,  Pa   and  Re   if  it  is  viewed  as  the  limit  of  a  continuous 
potential. 

To  illustrate  the  evaluation  of  integrals  of  the  type  given  by 
(4.16),  the  Boltzmann  term  will  be  considered.   In  order  to  emphasize  the 
difference  with  hard  spheres,  the  low  density  square  well  coefficients 
will  be  expressed  in  the  dimenslonless  form 


"o*  '  \k,c>r         xo  ■  \v,o/xoR  D£  =  \v,o/DoS         (4-17) 


where  the  bard  sphere  Boltzmann  transport  coefficients  are  given  by  (2.4), 
It  is  well  known  that  the  integrals  for  the  Poltzmann  coefficients  can  be 
reduced  to  exnressions  involving  two  dimensional  integrals  over  the 
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relative  momentum  and  impact  parameter.   Remembering  that  the  Sonine 
polynomial  approximation  in  the    lowest  order  term  in  equations  (A. 2), 
(4.8)  and  (4.12)  gives 


n*  -  an,n/no 


n         2    mB    aX,0M0 


r"D)£  =^aD,o/(nD)nS  •  (4-18) 


where 


2 

(/dp.f.(l)D..(l)Mn..(].))"1 


n,0         .2        J   '  l  nv    -    ij         i  ii 


2p 


,        =  — ffdpf(l)pMp)1.  (4.19) 

r,o       s      J     l   ov   "  ii   is1  li 


The  expressions  for  the  shear  viscosity  and  thermal  conductivit}'  reduce  to 
the  integral, 


n*  =  X*  =  [«C2'2)(T*)]  l  ,  (4. 20) 


while  the  self  diffusion  coefficient  can  he  expressed  as 


(nD)*  =  [nn'n(T*)]  '  .  (4.21) 
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kT 
Here   T*  -  — ,    while    the   J2-integrals    are   given    by 


';(T*)    =  --B-  C^)J/"   Jdp.dp.dr   Pl.[llm   9(  1  .?)eR(  1 2)f    ( 1  2)  ]p.  . 


h    lrTr1 

e+0 

(4.22) 


Q  (T*)   = ( — )  IHp   dp   dr  n    .  p    .  }      I  1 1 m   8(12) 

„ „      2     imr  J      1      2      l    h    lrTrL      „ 

3  2m  a  e+0 


x  ER(12)f0(12)(l+P12)]{pllP   } 


Tr 


(2  2) 

Taking  the  ft   '    integral  as  an  exanple,  it  is  reducer1  to  a  two 

dimensional  integral  through  the  steps  outlined  after  (4.14).   Using 
(4.16)  and  doing  the  integral  over  total  momentum  P  =  p   +  p   gives 


3  2a' 


Q*--'(T*)  =-— V^>  '  /dP  f0(p)  /dr^(p.r){plP.}Tr 


?  9  ■%, 

x  litr  {[r  o"ER(r,p)  2{p  p  }   ] 
e-t-0  r  Ra 


+  [o2e1'T*e^(r,p)2{p.pi}Tr]  + 


[rV  e1/T*eR(r,p)2{p.p.}T  ]    }.  (4.23) 

l  l  ri'TrJ„  -J 

Ra 


Thp  brackets  [•••]   indicate  the  functions  inside  are  to  ^e  evaluated  at 

it      +    +        _ 
initial  conditions  given  b}?  Irl  =  Ra  ,  a  ,  and  pa   respectively.   The 

n  •*-  ■+• 

operator  R(r,p)  is  the  resolvent  operator  for  relative  position  and 

/■      S/    n    f  2    9     a<Kr)  a   ^-1  * 

momentum  (i.e.,  R(r,p)  =  (e  H p.      -  -r-1 )    for  continuous 

m   i  3r .    3r .    3p . 
i      l      i 

potentials),  and  the  low  density  Histrihution  function  for  relative 

•       f    (    \         <%    v3/2      -Sp2/n 
momentum    is    f„(p)    =    ( — )  e 

0  Trm 
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To  evaluate  (4.23),  it  is  roted  that  for  initial  conditions  at  Ra 

-t- 
all  particles  scatter  free  so  that  p  <;oes  to  its  asymptotic  value  after  a 

finite  time;  therefore  the  relation  lim  eR(r,p)f(p)  =  11m  f(p;-t)  can  be 

used.   For  initial  conditions  at  a   and  Ra   for  scattering  states  the  same 

argument  holds,  however  there  is  a  cancellation  between  the  two  integrals 

so  that  the  sum  is  zero.   For  example,  for  initial  conditions  given  by 

Ra   on  the  forward  hemisphere  (p»o   >  ft)  which  miss  the  core  (the 

trajectory  leading  to  p   in  Appendix  C)  scatter  into  the  same  momentum 

state  for  particles  starting  on  the  back  hemisphere  (p_  with  impact 

parameter  greater  than  a")  and  thus  cancel  each  other  due  to  opposite 

signs.   Similarly  for  impact  parameters  less  than  a,  particles  starting  at 

Ra   are  cancelled  bv  particles  at  a     which  scatter  into  the  same  final 

state. 

It  has  long  been  argued  that  hound  states  give  no  contribution  to  the 

Roltzmann  coefficient;  for  the  square  well  this  can  he  explicitly 

20  +       - 

proven.  "      For  bound  trajectories  starting  at  a  and  Ra  ,  operator 

lim  eRjp.p.}    can  be  evaluated  by 
e+ft     1  1 

Vs   >  ^   -  +drfP"r){pJP1}Tr  llr1  ^{PiP^Tr 

Ra    ,a  e+0 

00 

>   1^1     _     +dr(pT){lim  e   /     dte~et{p  p   }      {p   (-t)p   (-t)}      . 

Ra   ,o  e+0  ft 


(4.74) 


Here  it  is  implicitly  assumed  that  C4.24')  is  to  be  integrated  only  over 

that  part  of  phase  space  giving  bound  trajectories.   By  explicitly  writing 

out  the  time  dependence  of  |p.p.r_   and  taking  the  Lanlace  transform   one 

1  l  i JTr 

gets  an  infinite  series  which  can  be  summed  into  a  closed  form  before  the 
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Unit  e  •*■  0  can  he  taken.   Taking  the  limit  then  proves  there  is  no 

contribution  from  bound  states. 

(2  2) 
Thus  the  ft  '        integral  is  given  hv 


0(2,2)       -BR   0   1/2  r,   *  ,-K  f     "  ,+  % 

fi     (T*)  =  — —  ( )  Idpf „(p)  J    dr  (p«r) 

„  2   mir  '    ( )    '  „  + 

Rm  P  a 


(4.25) 


x  lP,^}T>^}Tr 


where  p*  is  the  asvmptotic  value  of  p..   Since  the  dvnamics  is  fully 

1  3 

specified  hv  the  magnitude  of  p  and  relative  angle  (D*r),  (4.25)  can 
immediately  he  written  as 

fi(2'2)(T*)  =  --  (-^)3/2R2/   p7dn  fn(p)f    sin9  cos9  d0 
mm        '  _        0    '  „  + 
0  Pa 

,+•  ■*■    s2 
r(p*p*)  '    li 


P 


4       3^  (A'26) 


Changing  variables  to  reduced  momentum  y  =  /(3/m)  p  and  reduced 

impact  parameter  x  =  —  =  R  sin0,  and  denoting  p«p*/p   by  cosO(x ,  v") ,  gives 
a 

(?,  2) 

for  ft  the    result: 


00  2        R 

ft(2'?)(T*)    =   /   y7dy    e_V      /      xdx(l    -    cos20(y,y)).  (4.27) 


As  seen  in  Appendix  C,  if  x  is  less  than  a  critical  value,  a,  then 

cos0  is  given  by  the  function  L/(x,y);  for  larger  values  cosO  is  given  by 

(2  2) 
L9(x,y).   Thus  ft   '    may  be  written  as  a  sum  of  two  terms: 


(2.2) 
o/z>^(T*)  =  ^(T*)  +  T2(T*) 


with 
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Y        2  R 
I  (T*)  =  /   y7dye~y  /  xdx(l  -  L,  (x,y)) 

1        0         0 

»      _  2  a 
T  (T*)  =  /   V  dye  V  {/   xd-fl  -  T   (x,v)) 
T  0 

R  ? 

+  /  xdx(l  -  L2  (x,y))}.  (4.28) 


Fere  y    is  the  critical  momentum  where  all  trajectories  hit  the  core, 

_! , 

? 

y  T* 


2     1/2                1    1/2 
Y  =  (T*(R  -1))    ,  and  a  =  (I )    .   Similar  results  can  he  found 


for  ftn,1)(T*): 


n(1,1)(T*)  =  i1d(t*)  +  i2d(t*) 


Y        2   R 
I1n(T*)  =  2/   y5dye"y   /  xdx(l  -  L  (x,y)) 

o  o 


co  2       a 

i2d(t*)  =  2^  y_dye  y  M  xdx(1  "  Vx,y)) 

Y  0 


R 
+  /   xdx(l    -  L?(x,y))}.  (4.2") 


These  integrals  have  been  evaluated  over  a  range  of  temperatures  and  agree 

30  31 

with  the  results  of  Holleran  and  Fulhert   and  Brush  and  Lawrence. 

(II)       (2,2) 
Tahle  TV.l  gives  a  listing  of  the  ft   '    and  ft      integrals  for  R  = 

1.5.   Direct  calculation  of  (4.29)  for  high  temperatures  can  be  compared 

32  (1,1) 

to  perturbation  theory  for  the  square  well  gas    which  gives  the  ft 

integral  as  a  series  in  —  ,  (i.e.  ft   '    =  1  +  c  /T*  +  ...).   The  linear 
approximation  given  by  the  perturbation  appears  to  be  reasonable  only 
for  — r  <  .1:  at  — r  =  .1  the  error  between  ft   '    -  1  as  given  hv  (4.29) 
and  by  c  /T*  is  -  20%. 
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TABLE  IV. 1  COLLISION  INTEGRALS  FOR  THE  SQUARE  WELL  POTENTIAL 


T*  ST  '  ' (T*)  ft   '  ;(T*) 

.5  1.481  1.746 

1.0  1.331  1.477 

1.5  1.247  1.330 

2.0  1.192  1.247 

2.5  1.156  1.195 

3.0  1.132  1.159 

3.5  1.111  1.132 

4.0  1.096  1.113 

4.5  1.085  1.098 

5.0  1.076  1.085 

5.5  1.069  1.076 

6.0  1.063  1.066 

6.5  1.058  1.060 

7.0  1.054  1.054 

7.5  1.050  1.049 

8.0  3.046  1.045 

8.5  1.044  1.041 

9.0  1.041  1.037 

9.5  1.038  1.034 

10.0  1.036  1.032 

20.0  1.0174  1.0172 

50.0  1.0065  1.0061 

100.0  1.0031  1.0028 

200.0  1.00)5  1.0014 


CHAPTER  V 
EVALUATION  OF  THE  CROSS  TERMS  FOR  TFF  CLASSICAL  SQUARE  WELL  CAS 

The  first  density  correction  to  the  transport  coefficients  consists 

of  two  terms:   a  kinetic  piece,  A      ,  which  depends  unon  three-body 

kk ,  1 

dynamics,  and  the  potential  kinetic  cross  tern,  A     ,  which  depends  upon 

two-body  dvnamics.   The  calculation  of  A,  ,  ,  for  the  square  well  would  be 

kk,  1 

even  nore  difficult  than  for  hard  spheres.   However,  estimates  of  its 
value  can  be  found  from  molecular  dynamics  or  from  model  kinetic 
equations,  such  as  the  repeated-ring;  theory  for  self-diffusion  (see  for 
example  Figure  IT. 5).   The  cross  term  can  be  evaluated  in  detail 
for  n,  and  A,  allowing;  determination  of  the  importance  of  bound  states 
relative  to  scattering-  states,  as  follows. 


A.   Cross  Term  for  Shear  Vi.scositv 


From  equation  (4.6)  the  cross  term  for  the  shear  viscosity  can  be 
written  as 


n*   ,  =  lim   -   P  ,.  0,  Jdp1dp-dr(r.p.)6(12)£R(12) 


e+D  a  20mfi 

(5.1) 

x  fn(12)(l  -P12){P1JPli}Tr  . 

3  U.S. 
where  n*   ,  =  n  ,  , /o  Ti„    ,  and  the  first  Sonine  polynomial  approximation 
<|>k ,  1     <}>k ,  1     0 

has  been  used.   Proceeding  by  changing;  to  relative  and  center  of  momentum 
coordinates  leads  to 
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n*   ,  =  1.1  m  — ^-r^pr   /dp ,dn„d?(r  .p  .  W  1  2)ePd  2  ) 


e+0  a  20ni2 


x  yn)p-^(r)2{P.n.}Tr.  (5.2) 


Making  the  transforation  given  by  (A.  14)  gives 
_2  ... 


2,    ,  -3*(r) 

•    M. 

e+n  a  5m  fl 


"♦M  =  I!"  .3  J0(2,?.)  ^  f0(p)  /"r(p-D{[|r|  (r^e 

x  eR(r,p){p.p. }Tr]  + 


+  [|r|2(r.p  )e-3*(r)e^rlp){p.p  }Tr]   + 

Ra 


-  [|r|?(r.P.)e~34>(r)eR(r,p){p.P.}m  ]     .  (5.3) 

u      it  l  i  i'TrJ„  - 

Ra 


The  brackets  [...[    indicate  tbe  function  is  to  ^e    evaluated  for  initial 
r . 
i 

conditions  at  |r|  =  r. .   To  evaluate  tbe  operation  of  eR(r,p),  tbe 
integrals  are  separated  into  the  phase  space  for  scattering  states  and 
bound  states;  i.e., 


njk,i  -  "ik/8-8^  +  %w,i(b-"-)  (5-4) 


All  particles  wi  tb  initial  separations  |rl  =  Ra  will  contribute  to  the 

scattering  states;  for  initial  conditions  at     |rl  =  Ra   or  a   the 

■*■  "   2 

scattering  states  are  sjiven  by  tbe  condition >  e„.   Tbe  bound  state 

m       0 

contribution  is  given  by  complementary  portion  of  tbe  integrals  over 

[r|  =  Ra   and  a  .   For  tbe  scattering  states  eR(r ,p) {p. p .  }    in  tbe 

limit  e  ■*■  0  is  equal  to  tbe  asvmptotic  value  of  |p.  (-t)p  f-tH„  . 

1  i      l     JTr 
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Denoting  |p*p*}m  as 
1  1  i '  Tr 


llm  {P.("t)p  (-t)} 

t+oo 


iPifTr  ' 


P?P 


(5.5) 


the  scattering  states  part  can  be  written  as 


n*   .(S.S.)  =  — rr-rr  (— )   Jp  dp  e       Jd8sin8cos 


[r3[  r^  £|*    1  (r^  jj*^ 

p     p2     3   p   n2  P.a+ 


->-   ->--♦- 


+  el/T*[(r^*  ^    1  (r1p  p*  j 

L   p       2      3    p      2J  + 


-v-   ->--»- 


2 

2- 

P  a 

2 


^[(^(Jif)  -|(p)P^]   J, 
p  p  Ra 


(5.6) 


The  limits  of  the  integrals  for  various  initial  conditions  can  be 

identified  from  the  dynamics  given  in  Appendix  C,  and  the  functions 

*  ■*■       ■*--»- 
r  •  p5**      p  •  p^ 
( — ' — )  and    „      can  be  identified  with  one  of  the  K.  's  or  L. 's  there, 
p  2  11 

P 
Changing  to  reduced  variables  y  =  /g/m  p  and  x  =  Rsin0  allows 


(S.S.)  to  be  written  as  a  sum  of  six  integrals  corresponding  to 


<t>k,l. 
different  initial  conditions 


H*   /S.S.)  =    ^*    {J,  +  J_   +  e1/T*(J„  +  J.  +  J,  +  J,  )},     (5.7) 
4>k ,  1  (2,2)  L  In    2n  3n   4n    5n   6n 


where 


2   R 


6,  -y 
i   dye 

0  0 


r,   =  R  /   vdye  y   /  xdx(K.L.  +  -  K_) 

In         n  r>  *  *  3   0 
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00  ,    _  2   a  R 

J2n  =  R  /  y  dye  y  {/  xdxCK^  +  -  Kq)   +  J   xdxCK^  +  j  KQ)} 

Ry        _  ?  Ra 
J3n  "  '1/y^y  dye_y  '„  xdxCK6L6  "  K7L7  "  3  a  K0> 

J4n  =  1^*^   /nXdx(K6L6  -  K7L?  -  |  a2K<) 

RY    6     2  Ra 

%  =  "R  /i;  _y  ^e   Sn   xdx(KnLn  -  K8Lfi  -  3  a  Kfl) 

00  2   1 

^  =  "*  /1//^Hye"y  Un*d*(KllLll  "  K8L8  "  I  *2V 

Ra 
+  /  xdx(K9Lg  -  KgLfi  -  j  <x\)}.  C5.8) 


Here  y  and  a  have  the  same  meaning  as  in  Chapter  IV, 

y  =  (T*(R2-l))"l/2,  a  =  (1  +  l/(v2T*)1/2.   Also,  y  =  1//T*  is  the  value 
helow  which  all  trajectories  starting  inside  the  well  are  bound,  y  =  Ry 
is  the  value  where  all  trajectories  starting  inside  the  wel]  with 

impact  parameter  x  >  1  are  bound,  and  finally  x  =  Ra  (where 

2    ]  /2 
a  =  (]  -  l/(y  T*)    )  is  the  impact  parameter  where  for  a  given  value  of 

y  all  states  starting  inside  the  well  with  x  >  Ra  are  bound.   The  hound 

state  part  of  t\*  is  given  by  the  portion  of  phase  space  not  covered  by 

integrals  J_  -  J,  .   The  action  of  lim  eR(r,p)  for  bound  states  has  been 

evaluated  in  Appendix  C,  so  that  n*,  , (B.S.)  can  be  written 

<t>k,l 


,—  1  //¥*  2    1 

(*   c    ^        Mt  1/T*r//Vi      6,      -yV     .      1    ,,       2,1/2 

(B.S.)    =  ,         ,    e  {J  y    dye    J    J    xdx  -r   (1-x    ) 


♦k.r      ^  5J](2,2) 


+  /  _  y   dye   y   /     xdx  -kl-x    ) 

1//T*  Ra 
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l/T*  ,     2  1    ,       2  ,  .„ 

-R  J     y  dye  '  J  xdx  —  (1  -) 

n  0    J      RZ 

RY  ^  2  1  ,  2  ,  ,„ 
of  6j  ~y  f  j  1  r i  x  .1/2 
■R  J   y  dye    J   xdx  -r-  ( 1 -) 

1//T*  Ra  R 

RY  ,     2  R     ,       2 

-R  J   v  dye    J  xdx  -=■  ( 1 -) 

0  1     J      R 

-R  f  y%<r*'    J"  xdx  ^(l  -  ^x)1/2}.  (5.9) 

RY         Ra         R 


This  result  may  be  simplified  into  contribution  from  negative  energy 
states  and  positive  energy  states.   The  positive  energy  states  are  bound 
classically  but  not  quantum  mechanically.   Thus, 


R/tt       -1  /T* 
n*   .(B.S.)  =    n    7,    e  UL    (X.  +  X_  +  X,)  (5.10) 

*k>1         45J2(2»2)         !     2    3 


where 


v     t   r>3\  f      "  ,   -y  /,   2x3/2 

X    =  (-R  )  J    y  dye    (1-a  ) 

1T1  1//T* 

l//f*        2 
X2n  =  (1   }  ^     y  ye  Y 

0 

X,   =  /  Y  __  y6dye"Y  (l-R2aV/2.  (5.11) 

Jn     1//T* 


Here  X^  gives  the  contribution  from  negative  energy  states.   In  Table  V. 1 , 
(5.7)  and  (5.10)  are  given  for  several  R  values  and  various  values  of  T*. 
These  results  are  plotted  in  Figure  V.l  for  R  =  1.5  for  comparison  with 
hard  spheres  and  molecular  dynamics.   The  large  deviation  from  the  hard 
sphere  result  even  for  the  scattering  states  shows  an  important  dependence 


Figure  V.  1 .   Comparison  of  n*    and  molecular  dynamics   (AAA). 
Scattering  states  only  (-  -  ^TJ  scattering  states  plus  bound 
states  C ).   Horizontal  line  is  hard  sphere  limit. 
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on  potential  model.   Further,  there  is  a  significant  dependence  on  the 

hound  states  (i.e.,  approximately  307  at  T*  =  2.0;  using  the  Boltzmann 

o 
shear  viscosity  parameters  this  corresponds  to  -  300   K  for  Argon).   This 

important  effect  has  heen  seen  to  be  due  mainly  to  the  positive  energy 

bound  states. 

R.   Cross  Term  for  the  Thermal  Conductivity 

The  cross  term  for  the  thermal  conductivity  (equation  (4.8))  can  be 
written  as  the  sum  of  two  integrals: 


X\k,l  =I1  +V  (S-12) 


where 


3 
Ij  =  £im  —  g  (2  2)  /dPdpd?  <j>(r)  P±    eR(12)  fQ(12)  (s  (1 )  +  S^)] 
e+0  ma  30ft   ' 

3 

T2    =  ZilT1  — 3        (2   2)    Id*dVd*    CP'P)r.e(r,p)eR(12)f0(12)(s.(l)    +   S   (2)). 
e+0  mcf  30fi      ' 

(5.13) 


*  3  H .  S . 

Here  A  .  .  =  A.,  , /a  \n        ,  and  the  first  Sonine  polynomial  approximation 

$K,  1      <pK,  1      (' 

has  been  used.  The  second  integral  I2  can  be  handled  in  the  same  manner 
as  the  shear  viscosity  cross  term,  leading  to  equations  similar  to  (5.7) 
and  (5.10);  i.e., 


I2  =  I2(S.S.)  +  I2(B.S.).  (5.14) 


Here  the  scattering  states  contribution  is  given  by 
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VS'S->   ■   .*%.,!)    Kx  +  J2X  +  el/T   (J3X  +  J4X   +  J5X  +  J6X^ 


15  ST 


where 


Ju   =   R/JyV   e   y   /^xdx[2y2L4K4  +    (4y2  -   5)KQ] 


(5.15) 


.4,     -y'lf^-joJ 


J2X  =  RV  dy   e        {/0xdxt2v"L4K4  +   (4y"  "   3)KrJ 


+   /\dx[2y2L2K2   +    (4y2   -    5)K0]} 


J_.  =  /RY   y4dy  e  y  /Raxdx[2y2(L,K.  -  L,K_) 
JA     i//t*  U  h  h     /  / 


+  (4y"a-  -  5)K*] 


J4X   =  £/d*   e'y  /ixdx[2y2(L6K6   -  L^) 


+    (4y2a2   -    5)K«] 


JH  -  -<  _/*  «"y  70RV,[2y2anKu  -  ,,v 


+    (4y2a2   -   5)KQ] 


J6X  =   -<Yy4dv   e_y    {/Jxdx[2v2(TJuK11   -   L9K9) 


+(4y2a2   -   5)KQ] 
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+  /^xdx[2y2(L8K8  -  L9Kq)        +  (AyV  -  5)K  ]}.       (5.16) 


The  parameters  and  functions  in  (5.16)  are  the  same  as  those  in  (5.8), 
Similarly  the  hounri  states  contrihution  may  be  written 


VB-S->    =   1355^>2)el/T(XU  +  X2A  +  X3X)' 


where 


2 
X.,  =  (-R3)/°°  __y4dv  e~y  (v2-15/14)(l-a2)3/2 
iA         1//T* 

X2X  "  fl  ~  ^/("/^V^  e"Y  (y2  -  15/14) 

2 

X,,  =  JRy  _j4dy  e"y  (y2  -  15/14)(1  -  rV)3/2.  (5.18) 

JA     1//T* 


In  (5.18)  X«  gives  the  contribution  of  negative  energy  states  to  1^. 

The  first  integra]  in  (5.13)  is  of  a  different  type  from  that  found 
for  the  shear  viscosity,  and  it  is  identically  zero  for  hard  spheres. 
However,  since  the  operator  eR(12)  acts  on  a  function  of  the  energy  and 
momenta  only,  the  property   _  eR(12)g(12)  =     g(12;-t)  can  be  used  to 
evaluate  the  integral.   Doing  the  center  of  momentum  integral  gives 


2 

T  __3 fj^j"*"^    /    Nx/    \   -3$(r)   lira     3,        N/gp         3.  ,.    1Qs 

Tl    =  I    3n(2,2)  /^P^f0(p)*(r)e       e-0  eR(r'p)r^T-  ^      (5-19) 

The  combination  <()(r)  e       restricts  the  dr  integral  to  points  initially 
within  the  force  range.   The  action  of  the  resolvent  operator  on  the 
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2  2 

function  of  p   leaves  it  unchanged  for  bound  states  and  replaces  p   by 

2  2  r  + 

a  p  for  scattering  states.   The  dr  integral  is  seen  to  be  the  volume  of 

that  portion  of  the  well  where  scattering  occurs  (for  the  contribution  to 

Ij  from  scattering  states)  or  the  volume  of  that  portion  of  the  well  where 

the  initial  conditions  lpad  to  bound  states  (for  the  bound  state 

contribution  to  I,).   Changing  to  reduced  variables  x  and  y  shows  that  I, 

can  be  written  as 


\    --^^^-il^S.S.)    +   VB.S.))  (5.20) 


with 

2 
I    (S.S.)   =  /RT  _y2dv4  "   ay   >   e"y  F   (Ra) 
1//T* 

2 
+  ilyy2^    (f  "   <*V)    e"V   F2(Ra),  (5.21) 

and 

A    fl//T*   2  ,    .3  2,      -y2    .3 


IjCB.S.)    =  ^  j^11    y   dy(y  -   y^)   e   y      (RJ   -    1) 


+  iRT      y2dy(y-y2)  e^&R3  -  i)  -f.(r«)) 

2 

+  J^y^ycf  -  y  )  e_y  (y(*3  -  D  -  f2(RcO).         (5.22) 


In    (5.21)    and    (5.22)    the   functions    F,    and   F^    give    the   volume   of    the   space 
corresponding   to   scattering   for   the    square  well, 


F,(x)    ~(2  -   (2   +4>(1    -4>1/2) 
R  R 


.    -   2pM        x2    1/2  1  2.1/2 

+   2x   R(l -)  -   2x    ( 1    -  x    ) 

R 
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-}  (2  -  (2  +  x2)(l  -  x2)1/2) 

3  2       2 

F  (x)  =^-  (2  -  (2  +2-)(l  -^y)1/2) 
R       R 

2 
+  2x2R(l  -  ^r)1/2  -  \   .  (5.23) 

R 

* 
A  graph  of  X     is  shown  in  figure  (V.2)  for  R  =  1.5.   A  collection  of 
<j)k,l 

A*   ,  for  various  values  of  R  is  included  in  Table  V.2.   For  completeness 
<t>k,] 

a  table  of   n*   ^  and   A*   .  for  different  R  values  is  also  included 
k k , 0       kk , 0 

(Table  V.3).   It  should  be  noted  that  the  inclusion  of  the  bound  states 

leads  to  qualitatively  different  behavior  from  that  seen  for  the  shear 

viscosity;  this  should  be  contrasted  with  the  scaled  Enskog  and  modified 

Enskog  theories  of  chapter  II,  which  predict  the  same  behavior 

for  A     and  n.    . 
<pk,  1      (j)k,  1 


Figure  V.2.   Cross  term  for  the  thermal  conductivity   A*     ; 

(-  -  -)  scattering  states  only;  ( )  scattering  plus  bound 

states. 
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CHAPTER  VI 
DISCUSSION 

The  contribution  to  the  first  density  correction  to  transport 

coefficients  has  been  seen  to  consist  of  a  kinetic-potential  cross  term 

and  a  kinetic-kinetic  tern  involving  three  body  dynamics.   The  cross  term 

for  shear  viscosity  and  thermal  conductivity  have  been  calculated  in 

detail  for  the  square  well  potential  and  are  seen  to  have  a  strong 

temperature  dependence  both  in  the  scattering  states  and  the  bound  states, 

as  compared  to  the  hard  sphere  value  (Figure  V. 1  and  V.2).   The  kinetic 

terms  have  not  been  calculated  yet  for  an  potential  other  than  hard 

spheres.   However,  estimates  of  its  temperature  dependence  for  n* 

and  (nD)*  may  be  found  from  the  square  well  molecular  dynamics  data  of 

Michels  and  Trappeniers.  '     The  kinetic  part  for  shear  viscosity  ranges 

from  -  -1.1  at  T*  -    5.0  to  -  -1.8  at  T*  =  1.32,  compared  to  the  hard 

sphere  value  of  -1.26.   (The  uncertainty  in  n*   , f or  the  square  well  is 

large  enough  so  that  it  is  not  possible  to  say  whether  or  not  n*   ,  has  a 

maximum  value  and  approaches  the  hard  sphere  limit  from  above.)   The  self 

diffusion  value  is  seen  to  range  from  -  -1.5  at  T*  =  5.0  to  *  -5.0  at 

T*   =  1.5,  compared  to  a  value  of  -1.2  for  hard  spheres.   (It  is  seen  that 

the  first  density  correction  terms  and  the  Roltzmann  transport  coefficient 

all  have  similar  temperature  dependence.   This  may  account  for  the 

reasonable  agreement  found  between  scaled  bard  sphere  theories  and 

experiment.   )   As  already  noted  by  Micbels  and  Trappeniers,   the  sign  of 

the  first  density  correction  of  n  for  the  square  gas  will  change  sign  at 

around  T*  ^  3.0.   This  behavior  has  not  been  observed  in  experimental 
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results  on  simple  fluids.   However,  the  lack  of  data  in  the  low  density 

and  low  temperature  region  leaves  the  possihility  open.   For  example,  the 

33 
measurements  of  Haynes    for  the  first  density  correction  for  Argon  and 

Fluorine  above  the  critical  temperature  show  a  tendency  to  decrease  as  the 

temperature  is  lowered. 

As  mentioned  in  Chapter  II,  the  triple  collision  operator  as  defined 

by  (B.4)  has  been  calculated  in  detail  for  hard  spheres.   To  see  the 


connection  between  the  resul t  here  and  the  hard  sphere  calculations  of 


Hoegy  and  Sengers  ,  note  that  the  operator  M  can  be  written  (in  the 


classical  limit)  in  the  form 


M2y(l)  =  -|  f'Vl)  /dx2dx3{Rn1(123)(R(123)  -  RQ(123))W(123) 


I   R^CaXRCa)  -  R0(a))W(a)  -  I        I        Rn*  (a] )  (R^)  -  ^(c^)) 
a  a  a2*a\ 

-1 


x  W(0l)[fa  +  R0  (a2)(R(a2)  -  Rn(a2)W(a2) ] } 


x  f0(l)f0(2)f0(3)  (1  +  aP12  +  aP13)y(l).  (6.1) 


Here  a  denotes  the  pair  (i,i),  i,i  =  1,2,3  and  the  sum  is  over  distinct 

pairs;  f    is  the  Mayer  f-function.   For  the  special  case  of  hard  spheres, 
a2 

it  can  be  shown  that  the  operator  can  be  written  as  Mayer  f-functions 
times  operators  which  generate  sequences  of  binary  collisions. 


M0y(l)  =  -fn1(D  (dx0dxjy  f  f.T  +  I      f  T.R-T   +  ...}  (6-2) 

2        0     J   2   3lL  a  3  y   „    a  3  0  y      ' 

y        e,y 


where  the  T  's  and  T  's  are  the  binary  collision  operators  for  hard 
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spheres.   (See  ref.  (34)  for  a  discussion  of  the  properties  of  hard 

spheres  which  allow  a  simplification  of  this  result.)   In  this  way  the 

dominant  collision  sequences  can  he  identified.   For  hard  spheres  the  terra 

proportional  to  only  one  binary  collision  among  three  particles,  modified 

by  the  excluded,  volume  correction  due  to  the  Mayer  f-f unctions,  gives 

better  than  90%  of  the  triple  contribution  to  A     •   This  term  can  be 

identified  exactly  with  the  Enskog  correction  due  to  three  body 

interactions.   The  next  most  important  term  comes  from  two  binary 

collisions  among  three  particles  modified  by  the  excluded  volume  of  the 

non-interacting  particle.   Next  in  importance  are  the  'ring1  events, 

sequences  of  three  binary  collisions  among  three  particles,  while  all  four 

binary  collision  scattering  events  are  of  negligible  importance.   Due  to 

the  dominance  of  the  Enskoe  term,  attempts  have  been  made  to  derive  a 

similar  expression  from  (6.1)  by  a  binary  collision  expansion  of  the  full 

square  well  operator.   That  is,  an  expression  where  the  dynamics  of  the 

interacting  particles  is  modified  by  static  correlations.   However  the 

relationships  between  the  T-operators  and  the  Mayer  f-functions  for  hard 

spheres  which  allow  one  to  factor  the  f-functions  to  the  left  of  the 

binary  collision  operators  in  (6.2)  do  not  hold  true  for  other 

potentials.   In  equation  (6.2),  for  example,  the  first  term  is  obtained 

from  terms  like  f   (1  +  f  )f  t    ',    for  hard  spheres  T  f   =0,  but  not  for 
Y       Y   <*  3  Y  Y 

other  potentials.   The  Enskog  approach  is  also  seen  to  be  inadequate  for 

the  consideration  of  bound  state  effects,  as  seen  in  the  calculation  of 

29 
the  repeated  ring  theory  for  self-diffusion.     There  the  effects  due  to 

bound  states  are  a  major  contribution  to  the  first  density  correction. 

Thus  for  fluids  with  attractive  forces  it  may  not  be  reasonable  to  attempt 

to  separate  static  and  dynamic  correlations  between  the  particles.   From 
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the  above  discussion  it  is  seen  that  it  is  not  clear  how  to  derive  an 
Enskog-like  piece  from  the  triple  collision  operator  for  the  square  well 
potential  considered  here,  although  several  'Enskog-like '  models  have  been 
suggested  previously  for  evaluation  of  the  triple  collision  operator. 

The  dependence  of  the  self-diffusion  coefficient  for  the  square  well 
is  known  not  to  be  the  same  as  the  diffusion  coefficient  for  the  Lennard- 
Jones  or  for  experimental  data.   In  Figure  VI. 1  the  addition  of  the  square 
well  gives  poorer  agreement  with  experiment  than  the  Enskog  theory, 
although  the  decrease  in  the  self  diffusion  coefficient  for  the  square 
well  relative  to  hard  spheres  is  what  one  would  expect.   This  effect  and 
the  strong  temperature  dependence  found  for  the  cross  terms  may  he  due  to 
the  inappropriateness  of  the  square  well  as  a  model  for  real  fluids,  or 
due  to  using  classical  mechanics  to  calculate  the  bound  state 
contributions  to  the  transport  coefficients.   These  cases  will  be 
considered  below. 

For  the  first  case  it  has  already  been  noted  that  the  square  well 
parameters  for  a  fluid  as  determined  from,  P-V-T  data  and  those  determined 
from  transport  phenomena  are  sometimes  quite  different,  while  those  from  a 
smooth  potential  such  as  the  Lennard-Jones  model  are  close.   While  this 
indicates  that  one  may  not  be  able  to  map  fluid  behavior  for  the  square 
well  onto  experimental  results,  it  does  not  explain  why  there  would  be 
qualitative  difference  in  the  behavior  of  the  square  well  fluid  as 
compared  to  other  fluids  with  attractive  potentials.   High  temperature 
effects  due  to  the  hard  core  have  been  seen,  e.g.  the  second  virial 
coefficient  exhibits  a  maximum  for  real  fluids  due  to  penetration  of  the 
core  while  the  square  well  monotonically  approaches  its  hard  core  value. 
Another  effect  associated  with  the  discontinuous  potential  is  the 


Figure  VI. 1.   Comparison  of  nD  for  the  square  well  (AAA">,    Enskog 
theory  C ),  and  Lennard-Jones  and  experimental  data  ( ). 
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existence  of  high  energy  bound  states;  i.e.,  no  matter  what  the  angular 

momentum  between  the  two  particles  is,  there  is  always  a  region  of  phase 

snace  where  the  particles  are  bound  together.   It  is  also  well  known  that 

the  square  well  model  (nor  other  two  body  potentials)  cannot  match  both 

second  virial  and  third  virial  coefficients  of  experimental  data  due  to 

the  lack  of  three  body  forces.     In  fact  the  Lennard-Jones  and  square 

well  predictions  for  the  third  virial  coefficient  do  not  agree  (Figure 

(VI. 2)  even  though  the  second  viria]  coefficients  do.   However,  inclusion 

of  three-body  forces  may  not  resolve  the  discrepancy  between  the  square 

well  molecular  dynamics  and  experiment,  since  the  Lennard-Jones  molecular 

dynamics  data  for  self-diffusion  does  agree  with  experiment. 

The  close  agreement  between  the  results  presented  here  and  the 

molecular  dynamics  results  can  lead  one  to  argue  that  the  disagreement 

with  experiment  is  due  to  using  classical  mechanics  for  the  bound 

states.   Quantum  effects  would  be  most  notable  in  the  negative  energy 

states  which  are  hound  in  both  classical  and  quantum  systems,  and  in  the 

metastable  states  which  are  bound  classically  due  to  the  centrifugal 

barrier  but  not  bound  in  quantum  mechanics.   Figures  II. 2,  V.]  and  V.2 

both  exhibit  large  effects  due  to  bound  states;  for  tbe  shear 

viscositv,  n*   ,,  almost  all  the  bound  state  effects  are  due  to  the 
(f>k,l* 

metastable  states.   The  fraction  of  monomers  and  dimers  can  be  estimated 

35 
from  the  calculations  of  Stogryn  and  Hirschfelder'   (see  figure  VI. 3).   It 

is  seen  that  the  square  well  has  approximately  twice  the  fraction  of 

dimers  as  the  Lennard-Jones  for  the  same  state  conditions.   Furthermore, 

the  square  well  parameters  used  in  figure  (VI. 2),  found  from  PVT  data, 

have  a  much  shallower  well  depth  than  those  found  for  the  shear  viscositv 

(see  figures  1.1  and  1.2)  so  that  using  the  square  well  to  model  a 
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transport  property  leads  to  a  much  larger  fraction  of  ditners.   Unfor- 
tunately it  is  not  possible  to  predict  what  overall  effect  quantum 
mechanics  would  have  on  these  states,  so  the  use  of  quantum  mechanics  may 
have  no  effect  on  the  discrepancy.   (The  thermal  wavelength  X  is  typically 
of  the  order  .la,  so  one  would  not  expect  quantum  effects  to  he 
important. ) 

In  summary  this  work  has  presented  calculations  comparing  the  hard 
sphere  fluid  and  square  well  fluid  to  experimental  data  and  molecular 
dynamics.   It  is  seen  that  the  hard  sphere  fluids  cannot  descrihe  in 
detail  the  behavior  of  real  fluids  due  to  the  lack  of  an  attractive  well. 
The  success  of  the  repeated  ring  theory  for  self  diffusion  indicates  a 
generalization  of  it  for  a  more  general  type  of  transport  coefficient 
would  be  a  useful  approach  to  the  evaluation  of  higher  density  effects  in 
fluids.   (The  form  of  such  a  generalization  is  not  known  at  the  present 
time.)   The  square  well  calculations  agree  well  with  molecular  dynamics 
for  square  wells;  although  differences  still  exist  with  real  fluids.   The 
resolution  of  these  differences  would  imply  either  (1)  the  square  well 
fluid  behaves  in  a  manner  qualitatively  different  from  fluids  with 
continuous  potentials  or  (2)  that  it  is  necessary  to  use  quantum 
calculations  for  descriptions  of  the  behavior  of  dense  fluids.   These 
points  could  be  investigated  by  both  theoretical  calculation  and  low 
density  computer  simulation  of  continuous  potentials  for  comparison  with 
the  square  well  or  by  quantum  calculations  for  the  square  well.   Some 
investigation  of  molecular  dynamics  for  the  Lennard-Jones  model  has  been 
done.     The  slow  decay  of  the  time  correlation  functions  at  low  density 
appears  to  he  the  maior  obstacle  in  the  way  of  extensive  use  of  computer 
simulations  at  low  density,  the  latter  case  turns  out  to  agree  with  the 
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experimental  data,  it  would  indicate  that  the  square  well  potential 
emphasizes  the  difference  between  classical  and  quantum  systems  much  more 
than  continuous  potentials.   For  the  former  case,  it  would  indicate  that 
effects  important  here  are  sensitive  to  the  detailed  form  of  the 
potential.   It  is  indeed  possible,  that  both  effects  are  important,  and 
that  a  much  better  understanding  of  the  mechanisms  of  transport  phenomena 
in  dense  gases  is  necessary. 


APPENDIX  A 
CLUSTER  EXPANSIONS  OF  THE  TIME  CORRELATION  FUNCTIONS 

The  density  expansion  of  the  quantum  Green-Kubo  formulas  is  obtained 

1  2 
in  the  same  manner  as  for  classical  systems   : 

(1)  The  activity  expansion  for  n,  ty        ,  and  ty  is  obtained  from  a 
cluster  expansion  of  a  generating  function. 

(2)  The  expansion  for  n  is  inverted  to  give  the  activity  as  a  series 
in  n,  which  is  then  used  to  obtain  ty  and  \\i  or  a  series  in  the 
density. 

(3)  It  is  observed  that  the  reduced  distribution  functions  separate 

naturally  into  kinetic  and  potential  parts,  corresponding  to  J.  and  J 

and  it  is  seen  that  the  lowest  order  terms  in  \p   are  well  defined  from  the 

9 

cluster  expansion. 

(4)  The  kinetic  parts  are  seen  to  be  ill-defined  for  e  +  0  (or 

t  ■*■   °°),  hence  the  behavior  of  ik  is  to  be  determined  from  a  closure  the 

k 

BBGKY  hierarchy. 

(1)   Consider  a  time  correlation  function  of  the  form 


A  =  Wit.  ~   /dt  e  Et  <JAJA(-t)>  (A.l) 

e  ■*■  0 

v  +  °° 

with 


N 
Jh   =  (l  +  0     I      P   )  a(l)  +  |  I      b(ag).  (A.2) 

a* I  a*3 

"V 
and  J.(-t)  =  e     J.(0). 
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The  ensemble  average  is  given  by  (3.6)  and  the  Liouville  operator  is 

defined  by  (3.12).   P   is  the  permutation  operator  for  labels 
ag 

a  and  3,  and  a  =  0  for  self  correlations.  The  reduced  distribution 

(s) 

function  i\>  (l..s;e)  are  given  by  (3.9)  and  the  density  n  is  given  by 


p  =  v  +  m<n> 


£im  J_  y   Ne^  Tr 
N>0      1..N 


(A.  3) 


It  is  more  convenient  to  consider  the  cluster  expansion  of  the  reduced 
density  operators  defined  by 


N>s  *   B'--W  s+l..N 


rather  than  (A. 3)  and  (3.8)  separately.   Here  D  (l..N,a)  is  given  by 

-SH   aJ  C-t) 
DN(l..N,a)  E  N!  PN  X        e     e  (A. 5) 

N 
with  5(a)  =     I     ^y  D(l..N,a). 

N>0     1-N 

Here  a  is  a  c-number  parameter,  so  that  it  is  seen  that  Fv  J    is  related  tc 
\J;    and  n  by 


|(s)(1..8.t)  .aF(g)(l..B;,)| 

a=0 

(A.6) 

l 
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The  Husimi-UrseD  cluster  operators   ,  U(l.  .s  |r. .  t)  ,  are  the  most 
convenient  for  expansion  of  (A. 4).   They  are  defined  try 


n  (l..N;a)  =  U(l..s)  D„   (s+l..N;a) 

N  N-s 


+   I        U(l..s|i)  D     ({1  =  s+l..N|.1#i}  a) 
i=s+l  ' 

N 
+  I       U(l..s|ii)  DN_s_2  ({k  =  s+l..N|k#i,j}  a) 
i  <  i  s  1 

s  +  1 

+  ...   •  (A.7) 

By  setting  N=l,2,  etc.,  the  cluster  functions  can  he  found  in  terms 
of  D. *s: 


U(l)  =  Dj(l;a) 


D(l|2)  =  D2(12;a)  -  D1(l;a)D1 (2;a) 


(A. 8) 


Uf  12)  =  D  (12;a) 


U(12|3)  =  D3(123;a)  -  D  (12;a)  D  (3;a) 


Using  (A.7)  in  (A. 4)  gives 


nS  F(B>(l...;a)  -  J   (N.B^5(a)  {*(!..»)  Jr   D^Cj+1.  .N,a) 
N>s  s+1. .N 

N 

+  I       Tr  U(l..s|i)  Tr  n„   . f {i=s+l. .N|i*l},a) 
....  r . l   N-s-1 w  '       ' 

1=S+1  3  1 


=  eSUC]..s)  +  eS+1  Tr  U(l..s|s+1) 
s+1 

+  y  £S+?    Tr     U(l.. 8  I  8+1,8+2)  +  ... 

s+l,s+2 

00     £ 
..  nS  F(s)(l..s,a)  =  I      r~  Tr   U(  1 . . s | s+1 . .£)         (A. 9) 

£>s  U  S-       s+1. .A 

(S) 
From  (A. 9)  the  activity  expansion  for  n  and  ij>    can  be  found  from  (A. 8): 


s  t(s),,  v  v  £  n,  9U(1.  -s  fs+Jt)  I 

n    *     <1"8»t>-X   T^sTT   Jrp— W n 

£>s  s+1. . %  '  a=0 

n-vXin.?   ^    a?Tr    Tr    ufll2..*>|       •  (A-10) 

9U I 
Denoting  ■=■—     as  U',  it  will  be  assumed  that  the  cluster  operators  are 


9a 


a=0 


evaluated  at  a=0  from  here  on.   (Note  that  D  (l;a=0)  is  defined  such  that 
^  Tr  DjCl'.a-O)  =  1.) 

(2)   Inverting  the  equation  for  n  to  get  e  as  a  power  series  in  the 
density  gives 


o-  2 

e  =  n  "  V  +  -  T~   Tr   'TH|2)  +  ...  (A. 11) 

1,2 
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Substituting  into  the  activity  expansion  for  <jj    and  \\i  gives  (for 

convenience  the       will  be  assumed) 


$(1)(l;t)  =  U*(D  +  n  (Tr  U'(l|2)  -  TI'(1)  Tr   U(2|3)' 
2  2,3 


if>(2)(12;t)  =  U'(12)  +  n  (Tr  TT '  C  ]  2  I  3 )  -  2  u'^12)  Tr   U(3|4) 

3  3,4 


CA. 12) 


For  fluxes  of  the  form  given  by  (A. 2),  the  II'  and  IT  operators  in  (A.  12) 

are  (with  D(l..s)  =  D  (l..s;a)     ) 

S        'a=0 

Tl'(l)  =  D(l)  e"L(1)t  a(l) 


=  D(l)  a(l)  (A. 13) 


U'(l|2)  =  D(12)  e  L(12)t[(a(l)  +  oa(2))  +  ob(12) 
=  (D(l)  D(2)  a(2)  +  D(l)  a(l)  D(2)) 


U(l|2)  =  D(12)  e  L(12)t  -  d(1)  D(2) 


!J'(12)  =  D(12)  e  L(12)t[a(l)  +  aa(2)    +   ob(12)] 
(12|3)  =  D(123)  e"T'(l23)t[a(l)  +  o(a(2^  +  a(3) 


90 


+  o(b(12)  +  b( 13)  +  b(23))] 


-  D(12)  e  L(l2)t[a(l)  +  aa(2)  +  ob(l2)]  D(3) 


-  0(12)  D(3)  a(3) 


Therefore  the  first  terms  in  the  density  expansion  for  \J>    and  ty  are 

piven  by 


t(1)(l;t)  =  D(l)  (l  -  ^  Tr  U(2|3))  a(l) 

2,3 


+  n  Tr  [D(12)  e  L(l2)t  -  D(l)  D(2)]  (1  +  a?      )  a(l) 
2 


+  n  Tr  D(12)  e  L(12)t  0b(12) 


(A. 14) 


^2)(12;t)  =  D(12)  e"L(12)t  <*         2 


(l  -  ^  Tr  D(3|4))(l  +  a?      )  a(l) 

3,  h 


+  nTr  (D(123)  e^^K  _   ^    D(3)  e-L(12)t; 

3 


(1  +  a   P?  +  oP13)  a(l) 


+  D(12)  e  eL(12)t  (i  -  |  Tr   TJf3|4))  ob(12) 


3,4 


+  n  Tr  (D(123)  e  L(123)t  a(i  +  P^  +  P^)  b(12) 
3 


13    23' 


-  D(12)  D(3)  e  L(12;t)  abCl2) 
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The  integral  —  Tr   Ufa|3)  is  proportional  to  the  second 

a, 6 
coefficient : 


~   Tr   (D(12)  -  D(l)  D(2))  =  -  2nB?  (A. 15) 

1,2 


(3)   From  the  form  of  (A. 14)  the  potential  and  kinetic  parts  of  the 
reduced  distribution  functions  can  he  identified: 


^1}(l;t)  =  (1  +  2nB2)  D(l)  a(l) 


+  n  Tr  [D(12)  e  L(12)t  -  n(l)  D(2)]  (3+0  P   )  a(l) 


12' 

(A. 16) 


^2)(12;t)  =  D(12)  e  Ul2)t  (i  +  4nB2  )  (1  +  oP^)  a(  1 ) 


+  n  Tr  [DC123)  e"L(123)t  -  D(12)  D(3)  e-L(12)t] 
3 


x  (1  +  aP12  +  aP13)  aCl) 


and 


ijin)(l;t)  =  Tr  D(12)  e  L(12)t  0b(12) 

*  2 


^2)(12;t)  =  D(12)  e  L(12)t  (i  +  4nB  )  ab(12) 


r  ,    ,   -L(123)t  .    -    -  . 
+  n  Tr  [Dfl23)  e         (1  +  P  +  P   ) 

3 


-  P(12)  D(3)  e  L(12)t  a]    h(12).  fA.17) 
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Since  it  is  desired  to  work  with  the  Laplace  transformed  correlation 
functions,  the  transform  of  (A. 16)  and  (A. 17)  gives: 


^°(l;e)  =7  (1  +  2nB2)  D(l)  a(l) 


+  n  Tr  [D(12)  R(12;e)  -  -  D(l)  D(2)l  (1  +  oP,„)  a(l) 
„  e  J        12 


t(2)(12;e)  =  D(12)  R(12;e)  (1  +  4nB0)  (1  +  a?10)  a(l)         (A. 18) 
Tk  I  \l 


+  n  Tr  [(D(123)  R(123;e)  -  D(12)  D(3)  R(12;e); 
3 


(1  +  a?12  +  aP13)J  a(l) 


and 


4)^l)(l;e)  =  Tr  D(12)  R(12;e)  ob(12) 


i|/2)(12;e)  =  D(12)  R(12;e)  (1  +  4nB  )  ab(12) 


+  n  Tr  [n(123)  R(123;e)  (1  +  P  .  +  P..) 
3 


-  D(12)  n(3)  R(12;e)]  ob(12)  (A. 19) 

The  singular  nature  of  the  kinetic  parts  as  £  +  0  is  clearly 
indicated  in  (A. 18),  however  the  potential  terms  in  (A. 19)  (or  (A. 17))  are 
seen  to  have  well  defined  behavior  for  small  e  (or  t  +  °°).   This  is  due  to 
both  the  density  operators  D(ij..)  and  the  flux  terms  b(ii)  being  short 
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ranged  functions  of  position,  so  that  only  those  portions  of  phase  space 
where  the  particles  are  within  the  force  range  of  each  other  both  at  t  =  0 
(e  +  °°)  and  time  t  "•  —  will  contribute  to  the  integral  (A.l). 

(4)   From  Chapter  III,  it  is  seen  that  the  density  expansion  of  ij;,    and 

(2) 
ijj,    is  needed  for  evaluation  of  the  first  and  second  terms  of  the 
k 

density  series  for  A.   Since  (A.  18)  is  not  well  defined  for  e  ■*•  0,  the 
BBGKY  hierarchy  will  be  used  to  find  the  small  £  behavior  of  the 

correlation  function.   It  is  known  from  the  classical  results  that  the 

(2) 
explicit  singular  behavior  of  ik    at  e  +  0  can  be  removed  by  expressing 

(2)  ,  c    ,(D 

\li,    as  a  functional  of  ii,       : 
k  k 


^2)(12;e)  =  Y  ( 1 2 ;  e  |  ^  °  ( • ) )  -  (A. 20) 


Then  determining  i>~      or  a  series  in  the  density  from  the  first  BBGKY 
k 

equation  (Laplace  transformed)  - 


e   +  L(l))  ^(^e)  -  3£1)(1;t  =  0) 

(A. 21) 
(1). 


=  -  n/dx2  LT  ¥  [l2;e\ikU(')) 


-  will  give  the  kinetic  terms  in  a  series  well  defined  for  e  •*-  0. 

Equation  (A. 20)  easily  be  found  from  (A. 18)  by  solving  the  first  equation 

of  (A. 18)  for  a(l)  in  terms  of  ik;   (l;e),  and  then  substituting  this 

result  into  the  second  equation;  i.e.,  let  —  D(l)a(l)  =  I      a.n  ,  and 

E  i=0   n 

equate  powers  of  n  on  both  sides  of  (A. 18)  to  define  the  a. 's.   This  gives 


.(1)  =  (1  -  2nB?)  e  D  !(1)  ^l)(l;e) 
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-  neD  *(!)  Tr{[D(12)  eR,(12;e)  -  D( 1 )  D(2)] 


1  +  a?l2)   D  X ( 1 ) } ^^ U ( 1 ; e )  +  ...  .  (A.22) 


The  expression  of  i|>„    as  a  functional  of  ik    is  thus  given  by 


ikZ)(l2;e)   =  <P  (l2;e|^1}(.); 


=  D(12)ER(12;e)(l  +  a?     )  D  *(!)  ^l5(l;e) 


+  n  {D(12)6R(12;e)(l  +  aP   )  2B2  D  *(!)  ^l5(l;e) 


-  D(12)eR(12;e)(l  +  aP   )  D  l(l)  Tr  [d(13) 


eR2(13;e)  -  D(l)  0(3)]  (1  +  aP^)  D  *(!)  ^(l.e) 


+  Tr  [D(123)eR(123;e)  -  D(12)  D(3)  eR(12;e)] 
3 


(]  +  aP12  +  aP13)  D  l(l)   i|;[1)(];£) 


The  equation  above  can  be  written  so  that  the  terms  which  reduce  to  the 

classical  result  and  those  which  are  purely  quantum  mechanical  may  be 
identified : 

*(l2;e|^1)(»))  0(12)  eR(12;e)  (1  +  o?12)  D-1(l)  ^(^e) 
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+  n  Tr  [D(123)  eR(123;e)  -  D(12)  D(3)  eR(12;e)  (w(13)  eR(13) 
3 

+  W(23)  eR(23)  -  l)  (1  +  cr?12  +  oP13)]D^j)  tf'd.O 


+  n  D(12)  eR(12)  (1  +  aP^)  0(2)  D~*   t^  (l;e) 


+  ...   .  (A. 23) 


Here 


W(i.i)  =  D  l(i)   D  l(i)    D(1.1) 


-1 


(A. 24) 


ind   0(1)  =  2B?  +  D   (i)  Tr  (D(lj)  -  D(i)  D(j)). 

i 


From  (A. 24)  it  is  easily  seen  that  the  operator  0(1)  vanishes  for 
classical  mechanics.   The  first  two  terms  in  (A. 24)  give  the  classical 
result  if  the  replacements 


D(1..S)  +  fn(D  ..  fQ(S)  e 


,  I    4>(ii> 

i<j 


W(1..S)  ♦  e   Ki  (A- 25) 


R(l..s;e)  ■»•  (e  +  L  .  (l..s)) 
v     ex. 


,  +       +  ■*■       + 

and     Tr  +   Jdq.  dp.  ...  dq   dp 

'         1     1  S     £ 

1 .  .s 


Here  fn  is  the  Maxwell-Boltzmann  distribution,  <j>(ii)  is  the  interparticle 
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s   P„4     a       s 


potential,  and  L  ,   =   )   •  ■= +   }   LT(a3)  (note  LT(a3)  = 

cl.    **.   m    3q  .    50  I  I 

,.,  Qs   5  a=l         ai    a<3 


Mae,i   Hag,i 


for  continuous  potentials). 


Using  (A. 23)  in  (A. 21)  identifies  the  terms  in  Equations  (3.13)  to 
(3.19).   The  density  expansion  of  M  in  Equation  (3.15)  is  thus 


V(l)  =  e^O  D(!)  Tr  LT(12)  V12;e|D(1)  f(l)) 


M?f(1)  =  /^  D  l(l)    Tr  L^^^jelDd)  f(l)),  (A. 26) 


where  from  (A. 24)  ¥   and  1.   are  given  by 


f  (12;e|«)  =  D(12)  eR(12,e)  (1  +  aP n)   D  *(1)  (•) 


4-  (12;e|«)  =  Tr  [d(123)  eR(12^;e)  -  D(12)  D(3)  eR(12;e) 
3 


x  (W(13)  eR(13;e)  +  W(23)  eR(23;e)  -l) 


(1  +  o?n   +   0P13)]  D  L(l)  (•) 


+  J"1™      D(12)  eR(12;e)  (1  +  oP^)  0(2)  D  X(l)  (•)  (A. 27) 


As  mentioned  previously  the  operators  for  self  correlations  is  given  by 
a  +  0  above;  i.e. 


M,    f(l)   =     Mmn   D   !(1)    Tr   L_(12)    D(12)    eR(12;e)    f(l) 
is  e  ■+•  u  „       i 


M0    f(l)    =     *lnn   D   L(l)    Tr      {L_(12)    [d(123)    eR(123;e)    - 
2s  e  +  0  2)3lI 
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D(12)   D(3)    eR(12;e)    (w(13)    eR(13;e)   +  W(23)    eR(23;e)   -   l)] 


fC1)}    +      ^      D   1(1)    Tr   D(12)    eR(i2;e)    0(2)    f(l). 

e  +  0  2 


(A. 28) 


APPENDIX  B 
COLLISION  OPERATORS  FOR  THE  SOTTARE  WELL  POTENTIAL 

The  form  of  the  operators  given  in  Appendix  A  are  well  known  for 
continuous  potentials.   For  discontinuous  potentials  the  interacting  part 
of  the  Liouville  operator,  L  (i,i),  roust  he  defined  so  that  there  is  not 
explicit  reference  to  the  force,  or  the  collision  operators  must  he 
written  without  explicit  reference  to  the  L-j-'s. 

For  the  second  case  the  following  identities  are  used: 


s  -1 

I      L  (i,i)R(l..s)  =  -R  1(l..s)(R(l..s)  -  R  Cl..s)) 

Kj 

1 

Tr  L  (12)R(123)X(123)  =  Tr  L  (13)R(123)X(123) 

2,3  2,3 


Tr  L  (23)R(123)X(123)  =  0.  (B.l) 

2,3 


Here  the  e  dependence  of  the  operator  will  not  he  exhibited  explicitly. 
X(123)  in  (B.l)  is  any  function  invariant  under  interchange  of  particles 
and  R  =  (e+Ln)~   is  the  resolvent  operator  for  free  streaming.   From 
Appendix  A  the  first  and  second  terms  in  the  density  expansion  of  M(e) 


M.f(l)  =  0  V)  Tr  L  (12)eR(12)D(12)(l+aP  „)f(l) 

1  9 


M  f(l)  =  D  l(l)  Tr  {L  (12)[D(123)eR(123) 

2,3 
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-D(12)D(3)eR(12)[w(13)ERfl3)+W(23)eR(23)   -    l](l+aP     +oP     )]f(l)} 


+  D   l(l)    Tr   L   (12)D(12)eR(12)(l+oP      )Q(2)f(l).  (B.2) 

2 


In  the  limit  e  ■*■   0,  the  first  equation  in  (P. 2)  gives  the  Boltzmann 
collision  operator  for  quantum  systems  if  f(l)  is  diagonal  in  momentum 
representation.   The  second  equation,  for  M   ,  is  equivalent  to  the 
results  of  McLennan    for  the  triple  collision  term,  if  McLennan 's 
activity  expansion  is  converted  to  a  density  expansion.   Using  the  first 
equation  of  (B.l)  in  M   eives  the  following  result  for  M  : 


M  f(l)  =  -D  1(l)    Tr  R01(12)e(Rn2)  -  RQ(  12))D(  12)(  l+aP^f  (  1) .   (B.3) 


M  may  he  written  in  a  more  symmetric  form  hy  using  the  second  and  third 
identities  of  (B.l),  then  eliminating  the  L, 's  gives 


M  f(l)  =  -D  \l)    Tr  {j  R"1(123)e(R(123)  -  R  (123))D(123) 
2,3 


-RQ1(12)e(R(12)  -  R0(12)D(12)D(3)[WC13)eR(13) 


+  W(23)eR(23)  -  l](l+aP  2  +  oP23>f(l) 


-D  L(l)  Tr  Rn1(12)e(R(12)  -  R  ( 12) )D( 12) ( 1+oPj 2)0(2)f ( 1 )     (B.4) 
2 


The  interacting  part  of  the  Liouville  operator  Lj ,  is  defined  to  he 

the  commutator  with  potential,  d> .  . .   This  holds  true  even  for 

ij 

discontinuous  potentials;  however  one  must  he  careful  when  taking  the 
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commutator  with  momentum  operators  since  the  relation  |p   •]  +  —  - — 

f  L  a   J    1  3r 

a 

cannot  be  used.   For  the  classical  limit  this  implies  that  the  usual 
expression  for  LT, 


Li(l1,,-e     --Miii.' 


9r .  .     3p . . 

is  not  appropriate.   In  the  classical  limit  the  interacting  part 

of  the  Liouville  operator  L-j-(ii)  for  discontinuous  potentials 

generates  finite  changes  in  the  momentum  in  the  limit  t  +  0;  i.e. 

lim  e~      F(ii)  *   F(ij).   The  exact  form  for  the  interaction  term 

t+0+ 

depends  upon  whether  the  t  +■  0  limit  is  taken  from  the  positive  or 

negative  side,  whether  the  streaming  is  forward  or  backward  in  time,  and 

whether  the  probability  distribution  operator,  p,  are  acted  upon  by  the 

time  evolution  operator.   Since  this  work  does  not  explicitly  use  the 

interaction  operator,  it  will  not  be  considered  further. 


APPENDIX  C 
S0T1ARE  WELL  TRAJECTORIES 

The  dynamics  for  a  binary  collision  of  particles  interacting  via  a 

square  well  can  be  found  in  a  straight-forward  manner  using  conservation 

of  energy  and  conservation  ofargular  momentum  for  each  partial  collision 

at  the  core  or  the  outer  edge  of  the  well.   Depending  upon  the  initial 

conditions,  a  trajectory  must  be  one  of  seven  distinct  scattering 

trajectories  or  six  distinct  bound  state  trajectories  (Figures  C.l  and 

C.2).   In  Figures  C.l  and  C.2  and  the  following,  p  and  r  =  |r|r  denote  the 

+ 

initial  relative  momentum  and  position  of  the  particles,  and  p.  stand  for 

+ 
the  relative  momentum  of  the  particles  after  a  partial  collision  at  r. . 

r  i 

The  dynamics  are  to  be  calculated  backward  in  time,  and  most  functions 
will  be  expressed  in  reduced  units  of 


*   kBT 
temperature  T   =  


impact  parameter     x  =  b/a 


momentum  y  =  (—J    p  (C.l) 


where  b  is  the  initial  impact  parameter  and  -e  is  the  depth  of  the 
well.  To  solve  the  dynamics  it  is  useful  to  know  the  four  types  of 
partial  collisions  shown  in  Figure  C.3. 
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Figure  C.l.   Scattering  states  for  square  well  binary  collisions. 
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104 


105 


Figure  C.2.   Bound  state  trajectories  for  the  square  well. 
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II. 
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III. 
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VI. 


Figure  C.3.   Partial  Collisions  at  a    ,  Ra  and  Ro   for  the 
square  well. 


Ill 
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Using  the  conservation  principles  it  is  easily  seen  that 


(1)      For   p-r  >   0   and    r  ■*■  Ra 


HI)        +    .      a    A 

p  =     P    +    pgjT 


=  -ci  -  (~-/y/2  -  ci  -  &-)i/2 


=   (l   +    J/y2T*)1/2  (C.2a) 


(2)      For   p-r  >   0   and    r  +   a 


H2)        -  A 

p         =  p  -  pB  r 


2    1  /? 
$2  =   2    (1    -  x    )    '  (C.2b) 


(3)      For   p«r   <   0   and   r  +  Ro      (scattering) 


H3)        +  n,    A 

p  =   p  +   pB'    r 


x,2vl/2  r,         ,x    .2^1/2 


4  -  (i  -  (f>zrz  - « (i  -  <^n 


■  -a--H)1/2  (c-2c) 

y  t 
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(4)   For  p  •  r  <  0  and  r  +  Ro   (bound) 


+  (4)    +       * 
p    =  p  -  pg   r 


=  "  2  (1  -4)1/2  (C2d) 

R 


For  convenience,  the  scattering  states  and  bound  states  will  be  considered 
separately. 

1.   Scattering  states. 

Using  Equation  (C.2)  the  distances  a  particle  travels  between 
momentum  changes  are  easily  seen  to  be 

x  .2^1/2    r  x-il/2 


<,-*-<i/)"z-Ci-?) 


d2  -  R(l  -  <f/)>« 


a,  -  <>  -  &yn  -  (i  -  x2) 


\  -   "(1  -  (f)2)'72  (C3) 


Thus  using  Equations  (C.2)  and  (C.3)  one  can  find  the  eleven  position  and 
momentum  vectors  of  Figure  4  in  terms  of  the  original  vectors  p,  r: 


p  =  p  +  pg  r  |p  |  =  pa 


(C.4a) 


114 


P2   =   V   --S~J    P  +  P 


23, d        . 


'Pol     ■    P 


1  - 


Ra      ' 


-        2d2   + 
r    "   R^IP 


p3  -  (i  +  d^-)  p  +  pBj  [i  -j^  (r  --^; 


I P,      =  Pa 


(C.4b) 


=    (R- 


6ldl 


r p 

pa 


(C.4c) 


dlBl 


P4    =    t1+die2--ki(2   +   d132^    P 


dlBJ 


+  P   [2Bj  -  aB£  R  +  261B^1   -  -JLJ-  (2  +  d^')] 


IP/.  I    =  P 


(C.4d) 


t1  +  d,p2  -4r-^2  +  di^  ; 


-R^2   +   dl^P 


Note    that    in    (C.4c)    and    (C.4d),    3'    =    2    (l    --^r-)1^2. 


115 


+■  + 

Dr    =    P   -    P30    r 


Pel     =     P 


d    8*  0! 

P6  =  (l   -~f-)   p  +  p   (—  (1   +  d3B2) 


e2)  r 


i+  i 

I  Pa  I    =  Pa 


P?       =      [1       --R— J       P      +^ 


!p_      =   pa 


p     =  p  +   p$!    r 


pa 


2d.  3,' 
Pa    =   I ! 5— J    P 


!P9I    ■   Pa 


P10  =   ^   +   d33?J    P   "   pRB2   r 


P10l    =   P 


(C.4e) 


*  rl        d3*V       d3  + 

r6   =   ^R  +  T"J    r   "   RF  P 


(C.4f) 


1   *  3  + 

r7    =  R   r   "  R7  P 


(C.4g) 


(C.4h) 


3  + 
r.      =  R  r p 

10  p 


(C4i) 
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;11  =t1+  V2  --p^2  +  d3ei^  P 


-  P 


[R32  -  BJ  (l  +  d3B2)]  r 


11 


rll  '   Cl  +^32)  r  -_(2+d332);  (C.4i) 


In  the  calculation  of  the  transport  coefficients  the  functions 


■+-  4- 
P^ 

Li     2 
p 

(C.5) 
r*P,- 


i    P 

4- 

are  needed  for  the  asymptotic  states  p.,  1=2,  4,  6,  7,  8,  9  and  11. 
Defining 


«n=  (1  -X'J/- 


R2' 


%  =  (i  -  *Y'\ 


(C.6) 
2^1/2 


±     •+-        p*  r 

it  is  seen  that  for  initial  conditions  r  •+•  Ro   and  p«r  >  0,  =  K_, ,  for 

p     0 

■*-  •*■  +  * 

r  +  Ra  and  p«r  <  0,  ^—^-  =   -   Kn ,  for  r  +   a+   and  p.r  >  0,  ^^-   =  K'   and 
p        U  p      0 

■+■  " 

finally  for  r  -*-  o  and  p»r  <  0,  " —  =  -  K'. 

p       0 


117 


Therefore  the  functions  determining  the  Boltzmann  and  cross 
coefficients  for  scattering  states  are  given  by 


2d  ft  U, 

2d  R,  3,d? 

K2  =  ^1-^aLl)K0+2ei  Cl  "^  (C-7a) 


+  [231-Rag-  +  2R13^d1  -  —  (2  +  d^ )]  KQ 


,2 


d  g 
+  [2Bj  -  RaB'  +  2B1^d1  -  -^  (2  +  d^')]  (C.7b) 


d  b;         b; 

l6-(i  -rl^r^W-^^ 

K6  =  (l  -"-p-)    K'  +  (^  (1  +  d3R2)  -g2)  (C.7c) 


d  r;       b; 
d_g;  b; 

K7  =  -(l--fi)Ki+J-  (C.7d) 


S  -  !  "  *l    K0 


Kg  =  -  KQ  +  BJ  (C.7e) 
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2d  e; 

s  -  (•  — irj  -  % 

k9--(.  — ^)  +  «i 


CC.7f) 


l     =  [i  +  d3e2-^(2  +  d33;) 


-  [R32  -  &•  a  +  d3e2)]  kc 


Kll   =   []    +d3P2"4lf2  +  d3Pi)]   K0 


-  [RP2  -  BJ   (1  +  d332)]. 


(C.7r) 


Although  it  is  not  ohvious  from  the  form  given  here,  L,  ,  =  L^  and  L-,   =   Lq 

=  Lq.   The  L. 's  can  be  shown  after  some  tedious  algebra  to  be  equivalent 

to  the  functions  of  Marchetti  and  Dufty  (Equation  (B.31)  of  Reference  29), 

except  for  Ly,  L^  and  Lq,  where  the  first  term  in  their  expression  should 

2 
be—. 
R  a 

2.   Bound  states 

The  time  dependence  of  the  trajectories  shown  in  Figure  C.2  can  be 
expressed  as  an  infinite  sum  using  the  impulse  function,  I(t  ,t  ),  defined 
by 


Kt1,t2) 


(C.8) 
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For  the  trajectories  which  miss  the  core,  the  momentum  changes  are  spaced 
apart  by  a  time  T,  given  by 


T  =-i-f  (1-4)1/2  ,  (C9) 

4    /T*  Y      R2 


and  those  which  hit  the  core  by  a  time  To  given  by 


t  =-Ll  (r(i  -4)1/2  -  (i  -  x2)1/2).  (CIO) 

3    /T*  y        R2 


Note  that  the  definitions  of  the  other  reduced  variables  specifies  the 
2 

units  of  time  to  be  ( )    .   It  is  neither  convenient  nor  necessary  to 

£0 
find  the  momentum  or  position  in  any  time  interval  in  terms  of  the  initial 

p  and  r,  instead  the  functions  needed  for  the  bound  state  contribution  to 

the  transport  coefficients,  p*p(t)/p   and  r«p(t)/p,  can  be  found  directly 

from  the  initial  value  of  r«p/p  and  the  cumulative  angle  changes  in  p(t) 

due  to  collisions  at  the  core  and  the  outer  edge  of  the  well.   From  (C.2d) 

it  Is  seen  that  the  angular  change  in  the  direction  of  p(t)  for  collisions 

at  Ro  is  given  by 

+  .-*-        2 
cos  8   =  i      i+1  =  2~-  1,  (C.lla) 

P        R 

and  from  (C.2b),  collision  at  the  core  cause  an  angular  change  given  by 


cos  9„  =  2  x  -  I  .  (C.llb) 


Thus  the  time  dependence  of  r«p(-t)/p  and  p«p(-t)/p   for  the  six  classes 
of  bound  trajectories  can  be  written  in  the  following  manner 
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(I.) 


P'P 


2~  =      I      I(nT4,(n+l)T4)  cos  n  &2 


d     n=0 
r*Px 


~=  I     KnTA>(n+l)T4)  cos(n  92  -  y) 


,           , ,   x  v 1/2 
where  cos  y  =  O  ~) 

R 


(C.12) 


(II.) 


P*P 


II 


=  I     I(nT/it(n+l)TA)  cos((n+l)  %2) 


p      n=0 


r«p 


II 


=   £   I(nTA,(n+l)T4)  cos((n  +  1)  62  -  y') 

n=0 


u          i      /i    x  x  1/2 
where  cos  y  =  ~  (1 9) 

R 


(C.13) 


(III.)    P*PTII/P2  =  I   I(2nT  ,(2n+l)T  )  cos  (no) 
n=0 


+  I      I((2n-l)T3,2nT3)  cos  (na  -  8  ) 

n=l 


r»Pm/p  =  I      I(2nT3,(2n+l)T3)  cos  (no  -  y) 

n=0 


+  I      I((2n-l)T3,2nT3)  cos  (no  -  9?  -  y) 
n=l 


where  a  = 


2       r 


(C.14) 


(IV.)     p.p   /p2  =  I      I(2nT  ,(2ti+1)T  )  cos  (na  +  6.) 


IV 


n=0 


+  I      I((2n-l)T,,2nT,)  cos  (na) 
n=l 


3'   3' 


(V.) 
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r.p  00 

IV        V 

— - —  =  I      I(2nT3,(2p+l)T3)  cos  (na  +  6   -  y') 
n=0 


+  I      I((2n-l)T3,2nT3)  cos  (no  -  Y»),  (C.15) 

n=l 

■+■  -»- 

— Y  =      I      I(2nT3,(2n+l)T3)  cos  (na  -  8  ) 
p     n=0 


+  I      I((2n-l)T3,2nT  )  cos  (na) 
n=l 


r«p 

— —  =  I      T(2nT   (2n+l)T  )  cos  (na  -0-6) 

1      n=0  '  L 


+     I      I((2n-1)T  ,2nT  )  cos  (na  -  5) 

n=l 

2  1/2 
where  cos  5  =  (1  -  x  )  (C.16) 

->-  ■*■ 
P'PVT 
(VI.)     — ^  =  I      I(2n+1)T  )  cos  (na) 
p"     n=0 


+  I      I((2n-l)T3,2nT3)  cos  (na  -  6  ) 
n=l 


r'PV     v 

— —  =  I      I(2nT3,(2n+])T3)  cos  (na  -  6') 

n=0 


+  I      I((2n-1)T   2nT  )  cos  (na  -  0   -  6') 
n=l  J    J  1 


2  1/2 

where  cos  6'  =  -  (1  -  x")  (C.171 


From  Equations  (C.12)  -  (C.17)  the  time  dependence  of  any  function  of 

r«p(-t)  and  p»p(-t)  can  he  found;  the  I-functions  ensure  that  only  terms 
with  same  argument  in  I  are  non-zero-   The  series  can  he  transformed  term 


12; 


by  term,  and  the  limit  e  +   0  can  be  taken  after  the  series  is  summed  into 

a  closed  form. 

For  the  shear  viscosity,  the  e  +  0  limit  of  the  following  must  be 

found  for  the  bound  state  contribution  to  n  ,  , : 

<f>k,l 


£im     ,-«>  -etrf-r«p(-t)  yp«p(-t)-\    1  r«p- 

e  -»-  0 


r0  e-"[(£nizti)(P^ti)  _  i  i^p]  (c.18) 


Forming  the  product  in  the  brackets  and  taking  the  transform  gives  each 
case  (I  -  VI)  as  a  series  of  exponentials  and  cosines;  e.g.,  for  cases  I 
and  II: 

p.    °°        eT4    enT4 
(I.)     +  0  I      (l  -  e    )  e     cos  (n8  )  cos  (n6  -y)    -   —  cos  y 

n=0 

8  I     °°        "£TA   ~£nTA  1 

(II.)   eX,^0   H1  "  e    )e     cos((n+l)e2)cos((n+l)e2  -  y')  -  -^  cos  y' 
n=0 

(C.19) 

Similar  results  hold  for  cases  III  -  VI.   Although  it  appears  that  the 

series  are  zero  term  by  term  in  the  limit  e  +  0,  expressing  the  series  in 

a  closed  form  allows  one  to  identify  a  piece  which  is  non-zero  in  that 

limit.   The  series  can  be  summed  by  representing  the  cosines  as 

exponentials  and  summing  the  resultant  exponential  series.   By  example 

consider  case  (I): 


-eT.   °°   -neT      in6     -in6    i(n6  -y)  -i(nQ  -y) 

eliW0    {(1  -e   4)  I     e    4I(e   2+e    2)(e    2    +e     2    ] 
n=0 

1 
-  v  cos  y  = 
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tin.  Jfl  'e\    1   "  r  n(216  -eT  -iy)    -neT         -n(2i9  +£T  -iY) 
e  +  0{(l-e    J  4  Me  +  e     2  cos  y+e 

n=0 


cos  y  = 


-eT, 


,.    ,     4     -eT.+i(2e  -y)   ,       -eT. -i(26  -y)   , 
£im  rl-e     rr,     4     2    -\-l  .  fi     4     2  '  "\-l 


L  l  I    1 

+  —  cos  Yf  -  "o  cos  Y  = 

-eT,        -eT   cos(2e  -y) 
Aim  1-e     r   1+e  ,   1        1         ,  „„, 

e  ♦  0  2 I- ^T~ =2iTTJ  +  6  COS  T  =  6  C0S  Y    (C*20) 

l-2e     cos(26  -Y)+e 

Similar  results  can  be  found  for  the  other  cases,  giving  the  general 
result  for  the  bound  states  as 


tim    r°°  -£tr^r'p(-t)up'p(-t)\    I  r»p 
e  +  0   J0     Ll   p   n   p   >        3p- 


b   p 


This  surprisingly  simple  result  leads  to  the  bound  state  contribution  to 

ti  ,  ,  evaluated  in  Chapter  V. 
<j>k,l 
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